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Abstract 

A theoretical description of nucleon knockout reactions initiated by polarized electron scattering 
from polarized nuclei is presented. Explicit expressions for the complete set of reduced response 
functions (independent of the polarization angle) that can be experimentally obtained assuming 
plane waves for the electron are given in a general multipole expansion. The formalism is applied 
to the particular case of closed-shell-minus-one nuclei using two models for the ejected nucleon, 
including the final-state interaction phenomenologically with a complex optical potential and in 
the factorized plane-wave impulse approximation. Relativistic effects in the kinematics and in 
the electromagnetic current are incorporated throughout — specifically a new expansion of the 
electromagnetic current in powers only of the struck nucleon momentum is employed. Results are 
presented for the nucleus 39 K. 
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1 Introduction 



The completely unpolarized quasi-free (e, e'p) reaction has been systematically used to probe 
single-particle properties in complex nuclei, such as momentum distributions and spectro- 
scopic factors [1], and the reliability of the spectroscopic information deduced from these 
data is based on the weak dependence of the extracted momentum distributions upon the 
electron scattering kincmatical conditions of the experiment [2]. This is presumed to be 
especially true in the quasielastic region, namely at sufficiently high momentum transfer q 
and energy transfer lo <~ \/ q 2 + M 2 — M (with M the nucleon mass) where the process is 
dominated by the interaction between the electron and a single nucleon in the target, final- 
state interactions (FSI) are assumed to play as small a role as they can and the plane-wave 
impulse approximation (PWIA) should be expected to become at least roughly valid. Thus, 
to the extent that this simple picture applies, the knockout cross section can be factorized 
into the product of the electron-nucleus cross section and the spectral function, although in 
analyzing experiments one must usually take into account at least some aspects of the FSI. 

The last is certainly the situation when polarized electrons are used (but without hadronic 
polarizations — see later) , since the beam analyzing power involves a new observable — the 
fifth response function — that vanishes in PWIA and accordingly depends fundamentally 
upon the presence of FSI [3, 4]. Something similar happens when polarized nuclei are used 
as targets. In that case a total of nine classes of structure functions can be obtained from 
(e,e'N) reactions, specific ones of these being zero in PWIA [5]. This means that some of 
the polarization observables are expected to be sensitive to specific aspects of the reaction 
mechanism that are less pronounced in the unpolarized cross section. It is then of interest 
to study which are the relevant observables in these kinds of reactions and to explore their 
importance and sensitivity to different aspects of the modeling in the kinematical ranges of 
interest. 

In this paper we focus our attention on exclusive quasielastic scattering of polarized 
electrons from polarized nuclei. First, based on the general formalism for (polarized) coin- 
cidence reactions presented in Ref. [4] we specialize to the reaction A(e,e'N)B and, as an 
extension of our recent work on the inclusive reaction [6], we introduce in Sec. 2 a set of 
reduced response functions which characterize the polarized cross section and which do not 
depend on the polarization angles, but only on the momentum and energy transfer (q, lo) 
and the missing momentum p of the nucleon, having fixed the missing energy to knockout 
of a nucleon from a specific shell. To do this, we perform a multipole expansion of the 
nuclear states and electromagnetic current. This involves a sum over third components of 
the angular momentum that is performed analytically. The resulting exclusive cross section 
and hadronic structure functions also depend on the polarization and emission angles, which 
we can write explicitly as a linear combination of spherical harmonics. In particular, a mul- 
tipole expansion in terms of spherical harmonics of the polarization angles is obtained. The 
coefficients in this expansion are directly related to what we call angular reduced response 
functions. These functions constitute the basic ingredients which enter in the total cross 
section. We illustrate in Appendix A how to obtain the final multipole expansion for the 
particular case of the longitudinal response. 

Although we shall see that the number of angular reduced response functions is very 
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large, too much so to hope for experimental extraction of all of them in the general case, it 
is of theoretical interest to study their sensitivity to specific details of the reaction mechanism 
with the goal of deciding which are the most relevant to measure. This is the subject of the 
present paper. 

We explore these ideas within the context of the shell model. In Sec. 3 we specialize the 
general expressions to the particular case of nucleon knockout from closed-shcll-minus-one 
nuclei, leaving the daughter nucleus in a discrete state which is described by two holes in 
the core. We apply the formalism to the 39 K nucleus which is described as a d 3 / 2 hole in 
the 40 Ca closed shell core. Our theoretical framework addresses the final-state interaction 
of the ejected nucleon with the residual nucleus by solving the Schrodingcr equation with 
a phenomenological (complex) optical potential. The main result of this section is that in 
the extreme shell model the reduced nuclear responses are proportional to the polarized 
response functions for a single polarized particle in a shell. This is more or less obvious for 
one-particle nuclei, but not for the case of one-hole nuclei, because the involved shell and 
target nucleus can have different angular momenta, and the final daughter nucleus remains 
in a two-hole state with definite total angular momentum and then both holes are partially 
polarized. We prove the above result analytically for hole-nuclei by performing the sum over 
the total angular momenta of the final hadronic state composed of the daughter nucleus plus 
the ejected particle (see Appendix B). 

One of the aims of our systematic investigation of spin observables initiated in this paper 
is to determine whether momentum distributions of polarized nuclei (in particular nuclei 
near closed shells) may look different from those in the unpolarized case. In fact, single- 
nucleon knockout from polarized nuclei can provide a new probe of the spin-dependent 
nuclear spectral function which represents the probability of finding a nucleon in the target 
with given energy, momentum and spin projection. Then, spin observables can be used 
to explore the spin distribution of the single-particle orbitals, from which one can extract 
valuable information about the complete spatial distribution of the orbits. The general 
framework for such studies in PWIA was provided in Refs. [5], [7], and it was applied in 
Ref. [6] to the study of the inclusive polarized responses of one-hole nuclei. Accordingly in 
Sec. 4 we introduce the formalism needed for the particular case of one-hole nuclei. 

The factorized PWIA gives us a very clear picture of the initial state physics, but the 
final-state propagation of the outgoing nucleon has to be taken into account. In the present 
work we are able to go beyond PWIA and to determine to what degree the FSI effects can 
obscure the extraction of those single-particle properties. An exploratory study to set the 
scale of possible measurable quantities including the FSI was performed in Ref. [8] . 

It is important to note that the quasi-free conditions that favour the validity of the 
PWIA require a high value of the momentum transfer q which prohibits the usual non- 
relativistic expansions of the electromagnetic current in powers of q/M. Furthermore, for 
these kinematics the ejected nucleon becomes relativistic. Thus, in this work we use a new 
approximation to the on-shell relativistic one-body current that was tested in Ref. [9]; this 
involves an expansion only in powers of rj = p/M, but not in k = q/2M or A = ujjlM. In 
addition we incorporate relativistic kinematics throughout in the calculations. 
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2 General formalism for exclusive electron scattering 



2.1 Cross section 

We begin this section by considering the scattering of polarized electrons from polarized 
targets. The general formalism for coincidence electron scattering in which all of the particles 
involved can be polarized has been presented in Rcf. [4]. In the present treatment we have 
re-worked this formalism keeping in mind the specific process in which only the initial 
particles are polarized and written down the general equations in a way that is oriented 
towards studies of nucleon knockout reactions. For example, for the shell model applications 
considered here, it is convenient to use jj coupling instead of LS coupling as was used in 
Ref. [4]. Also, following the scheme introduced in Ref. [10], we make an expansion of the 
nuclear structure functions in a basis of spherical harmonics of the polarization angles, in 
that way introducing a set of reduced response functions that completely determine the 
process. 

Insofar as the electron scattering process is concerned, as discussed in Ref. [4] we limit 
our attention to the one-photon-exchange or plane- wave Born approximation (PWBA). The 
four momenta of the incident and scattered electrons are labelled = (e e , k e ) and K'^ = 
(e e , kg), respectively. The four-momentum transfer is given by = (u, q). We work in the 
laboratory system, where the initial nucleus is at rest in the state denoted by \A). After the 
interaction a proton or neutron with momentum direction p' and energy E' is detected and 
the (undetected) daughter nucleus is left in a discrete state \B). We can assume that when 
the nucleon is detected in coincidence with the scattered electron, it is essentially at infinity, 
and therefore on-shell, so that the magnitude of the asymptotic value of the momentum p' 
is related to the energy by E' = y V 2 + M 2 , as usual for free particles, although at short 
distances when the effects of the FSI arc important, the momentum of the ejected nucleon 
is not well-defined. Therefore, we define the momentum p' of the ejectile as the well-defined 
asymptotic value and so can label the nucleon wave function as |p') = \E',p'). 

On the other hand, later in using the formalism in the shell model where recoil is not 
handled properly we make no attempt to extract the center of mass motion from the nuclear 
wave function. Note that this is not a real limitation of the general formalism that follows, 
although a proper treatment of the problem is far from trivial. Indeed, an artificial extraction 
by hand of the center of mass motion in the shell model is known to generate spurious effects 
that only can be controlled in very specific cases, in particular not in the conventional 
continuum shell model. Our expectations here in applying the shell model are that the 
recoil effects in the quasi-free regime are not especially important for nuclei as heavy as 39 K. 

Also, for the moment only the electron is treated relativistically. Later we will approxi- 
mately include the relativistic effects in the ejectile nucleon by using relativistic kinematics 
and an appropriate electromagnetic current operator. Therefore, the S'-matrix element for 
the above process can be written as 

S fi = -2m5{E' + E B -E A - u)-^-(K'h'\j^0)\Kh)(p'm s , B\ J"(q)|A), (1) 

where h and h' are the helicities of the (ultra-relativistic) initial and final electrons, Ea and 
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Eb are the energies of the initial and daughter nuclei, Q 2 = u 2 — q 2 < is the spacelike 
four- momentum transfer, a is the fine structure constant, j^(0) is the electromagnetic cur- 
rent operator for the electron, m s is the spin projection of the ejectile, and J M (q) is the 
Fourier transform of the nuclear electromagnetic current operator. The differential cross 
section in the laboratory system is proportional to \Sfi\ 2 and to the electron- nucleon phase 

space ^-pd-E'dp', where dp' = dip' d cos 8' is the solid angle differential corresponding 
to the emission angles (6',(f>'). Performing the sum over final undetected electron hclic- 
ity and integrating over the nucleon energy E' using the energy conservation condition 
E' = Ea + u> — Eb, the resulting cross section is proportional to the contraction r q llI/ W lxv of 
the leptonic and hadronic tensors rj^ and W^ u given in Ref. [10] and below, respectively. 
In practice, the ejectile energy E' is detected, so the integration over the energy corresponds 
to experimentally suming up all of the events that contribute to the selected peak in the 
missing energy spectrum. 

Finally, the contraction rj^ v W^ v is performed as usual in the extreme relativistic limit 
for the electron and in the coordinate system with the z-axis along q and the x-axis in the 
scattering plane. Using current conservation, the result can be written [10] 

[ v L K L + v T n T + v TL K TL + v TT K TT + h{v T ,Tl T ' + v TL ,Tl TV ) (2) 
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= E + /iA, (3) 

where om is the Mott cross section, The electron kinematic factors vk are defined in Ref. 
[10]. The six exclusive nuclear response functions TZ K are the components given in Eqs. (9- 
14) below of the hadronic tensor for the exclusive process defined by 

W» v = <p , m a B|J"(q)|^)*<p , m a B|J"(q)|A) ) (4) 

m s M B 

where we sum over the undetected polarizations of the final particles. For a treatment of 
the general case where all of the particles are polarized, see Ref. [4]. Here we have followed 
a somewhat different derivation that is more appropriate for the particular application to 
the shell model. 



2.2 Polarized response functions 
2.2.1 Hadronic final states 

Now some more detailed specifications of the nuclear states themselves are in order. The 
initial nuclear state is assumed to be 100% polarized relative to an arbitrary axis of quan- 
tization whose orientation is specified by the angles f2* = (9* , <p* ) relative to the coordinate 
system defined above. We also denote by Jl* the unit vector pointing in the f2* direction. 
Therefore, 9* is the angle between q and f2*; note that the conventions here have <j)* with 
respect to the electron plane, whereas Ref. [4] starts with reference to the hadron plane. 
We express the nuclear wave function in terms of state vectors defined with respect to the 
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z-axis: 

\A) = \JiJi(sr)) = R(4>*,e*,o)\JiJi) = J2v<^ h (n*)\J*M t ). (5) 

Mi 

Note also that the three Euler angles (</>*, 9*, 0) of the rotation matrix are different from 
those specified in Ref. [4]. However, since any rotation around the SI* axis does not modify 
the polarization vector, the resulting observables are exactly the same in the two approaches. 

The final hadronic state corresponds asymptotically to a nucleon with energy E', mo- 
mentum direction p' with angles (9', </>') referred to the same coordinate system, and with 
(unobserved) spin projection m s , together with an (undetected) residual nucleus in the dis- 
crete state \B) (that in the general case is recoiling with momentum P B ). In analogy with 
the plane-wave expansion in spherical harmonics, we assume that a similar expansion can 
be performed for the hadronic states 

|p'm a B)= J2 i%*M{v')(\™slM\ 3 m){3mJ B M B \J } M } )\{E'l 3 )J B -J } M } ). (6) 

IMjmJfMf 

Here \(E'lj)J B ; JfMf) is a continuum nuclear state having asymptotic behaviour 

\(E'lj)J B ;J f M f ) \a Enj \J B )\ , (7) 

L J J f M f 

where a En - creates a single nucleon with energy E' and angular momenta (Ij). Note that 
the condition in Eq. (7) only has to be verified asymptotically. Therefore we do not make 
any assumption about the particular nuclear model to be used until later, keeping the 
following steps in this section completely general. However, it may be helpful to anticipate 
the modeling used later by noting that any particular theory of the reaction mechanism 
has to be able to provide a suitable form for the continuum states \(E'lj)J B ; JfMf) . In 
particular, in the simple case of the extreme shell model, Eq. (7) is assumed to be true not 
only asymptotically, but also at short distances. In the simplest case of the PWIA the radial 
wave function is proportional to a spherical Bessel function and then Eq. (6) reduces to a 
plane wave times the daughter state \B). As for the normalization of the states, it is 

((E'l'f)J' B ; J' f M' f \(Elj)J B ; J f M f ) = 6(E - E , )6 u .5 Jj .6j B j. B 8j fJ . / 6 M/M . . (8) 

This relation, together with the asymptotic condition in Eq. (7) uniquely determines the 
nuclear states. 

2.2.2 Multipole analysis 

In order to perform a multipole expansion of the responses, it is convenient to write them 
in terms of the spherical tensor components of the current [4] : 

TZ L = W 00 = J2p*P (9) 
n r = ^ XX + W W =^[|J_i| 2 + |J+i| 2 ] (10) 
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U TT = w yy_ wXX ^^[J^J+l + J* +1 J-l] (11) 

TZ TL = V2(W 0x + W x0 ) = - ^2Re[p* (J +1 -J_l)] (12) 
ft Ti ' = iV2(^ 0j/ -VK J ' ) = -^2Rc[p*(J +1 + J_ 1 )] (13) 
ft T ' = J (^- VK ya; ) =^[|J+i| 2 - |J_i| 2 ]. (14) 

Now we perform the usual multipole expansion for the charge and transverse current oper- 
ators: 

p(q) - V4^^ i J [J]M J0 (g) (15) 
J 

J m (q) - -a/2^]Tz j [J] 

J 

where to is the index of the spherical components and Mj m , Tj l m and T™^ 9 are the usual 
Coulomb (CJ), electric (EJ) and magnetic (MJ) multipolcs. Throughout we use the notation 
[J] = V2 J + 1 for any angular momentum variable J. 

Using the above expansion, the nuclear tensor can be written as a sum of factors of the 
type 

B$r= J2 (p'm a ,J B M B \f'j, ml \A)*(p'm a ,J B M B \f Jm \A), (17) 

m s M B 

where T and T' are any of the C, E or M operators. For instance, in Appendix A we show 
the expression for the L response given in terms of the B°p \j coefficients (see Eq. (132)). 
The parameters to or to' take on the value when the electromagnetic operator involved 
in the corresponding response is the charge one, while they take on the values 1 or -1 when 
transverse current components occur. Inserting the expressions Eqs. (5,6) of the initial and 
final nuclear states, respectively, into Eq. (17), we need to evaluate the following multiple 
sum 

m 3 M B I'M' j'm' p J' f M' f lMjm p JfM f Mi Mi 

x<^m s rM / |j / m;)0X J B M Bl4 M /><5 m - JM l^>0'^^M B |J / M / ) 
x ((l'j')JB , JfM' f \f' J>ml I JiM!)*{(lj)J B , JfM f \f Jm | J,M 4 ) . (18) 
The procedure to simplify the above expression is then the following. 

1. Write the product of the two spherical harmonics YiuYy M , as a unear combination of 
spherical harmonics Yji m'{v')- Therefore the index J' below refers to the multipole 
expansion in emission angles. 

2. Rewrite the dependence on the polarization direction in the product of two rotation 
matrices V^,j {VL*)V^'j,{VL*) by developing this in spherical harmonics of the po- 
larization angles Yj_m(Q*), and where the coefficients in the expansion also contain 



Tf m (q)+mfZ 9 (q) 



(16) 
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as a factor the Fano tensor for 100% polarization given by fj = (JiJiJi — Ji\JQ) 
(see, for example, Ref. [4]). Therefore, the indices J M. in the equations and reduced 
response functions defined below refer to this multipole expansion of the polarization 
dependence. 

3. Use the Wigner-Eckart theorem in the matrix elements. 

After these steps, all of the dependence upon third components is via 3-j coefficients, 
and the corresponding sums can be performed analytically using Racah algebra, but here 
we only quote the final result. Also, following Ref. [4], in order to simplify the notation it is 
convenient to define multiple indices cr, a' by a = {l,j, Jf, J} and a' = {l',f, J'f, J'}, and a 
coefficient 

$*>a(J, J',L; Jt -» J B ) 

As the values of Ji and Jb are fixed, we shall usually not write them in the argument of 
"JVcr- We quote here an important symmetry property of this coefficient under exchange of 
the two indices a 1 a that can be easily verified: 

*MJJ'L) = (-l) J+J ' +L $a*<(JJ'L). (20) 
With these definitions, the general expression for £?™j™ becomes the following: 

mWr = E E EE i ''" lp w'+j-(- 1 ) m /j^( fi, ) y j'(p')W 

IjJf l'j'J' f J J' LM 

x*mjj'l)(£, : J m m ) (nfy\\jy(f\\tj\w, (2i) 

where we define parity projectors Pf = (1 ± (— l) J )/2 and use shorthand for the generic 
final states, |/> = \{E'lj)J B ;J f ) and \f) = \{E'l'j')J B ;J' f ). 

Now, in order to write compact expressions for the different response functions it is 
convenient to define the following Coulomb, electric and magnetic multipole matrix elements, 

C a ee (f\\Mj\\Ji) (22) 
E a ee (f\\ffU) (23) 
M a ee (f\\ifr a \\Ji), (24) 

that in general are complex functions only of (q, u>), since the asymptotic value of the energy 
is given by E' = Ea + oj — Eb, and both Ea and Eb are fixed. The response functions 
are defined through linear combinations of the real and imaginary parts of the following 



[J] [A [j] if] [J f ] [J' f ] m [L] ^y+JB+Jf+^+j+j 

J J' L 



J J' \ f J J' 

1 n 



2 



J, J'f Jfl J I * 'ji 3 3 , <•») 
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quadratic products of the multipole matrix elements: 

(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 

so it is convenient to 

write the angular dependence in Eq. (21) as a sum of real plus imaginary parts. Thus we 
define the real functions Ajj>lm{Q*,p') and Bjj'lm(Q*,p') as 

Ajj, LM (sr,p') + iBjj, LM {sr,p') = [Yj(n*)(g>Yj,(p')} LM . (32) 

Finally we also define two parity functions 

&j = (-1) {J '- J)/2 P} +J , Oj = {-ir- J+1),2 Pj, +J . (33) 

After these preliminary definitions and following a procedure similar to that in Ref. [4] 
(summarized in Appendix A for the particular case of the L- response), we are then in a 
position to write the exclusive responses in the following relatively compact way: 





— 






= K'Ea- 


f M* a ,M a 




= E*,M a 


- M*,E a 


nTLl , ■ tTL 1 


= c;,E a 




D TL2 | ■ T TL2 
K'a + *V<r 


= c;,M a 




pTTl , • tTTI 


= E* a ,E a - 


- M*,M a 




= E*,M a 


+ M*,E tr , 


and Ie> a are real functions only of (q,ui) 



1l L = ^Y,fj{ P J A ^^^ W JJ'L-PjBjj'wW^-) L } (34) 

JJ'L 

K T = 4tt£ fjP^ +L {P^Ajj, L oW^ L -PjBjj ILO W^-l} (35) 

JJ'L 

TZ T ' = AnJ2fj P J' + L{PjAjJ'LoW^j7l-P+Bjj, L oW^jtl} (36) 

JJ'L 

U TL = A-k ^2 fj{PjAjj'LiWjjty—PjBjj> L1 Wjj^} (37) 

JJ'L 

U TL ' = ^ J2 fj{ P J A JJ'^ W JJ'L ] -PZBjj, L1 W%P} (38) 

JJ'L 

U TT = 4nY / fj{PjAjJ'L2W^, + L ) -P J Bj J/L2 W^- L ) }. (39) 

JJ'L 

The reduced response functions Wjj^ l (q, u>) contain all of the information about the nuclear 
structure and reaction mechanism and they are the following real quadratic forms involving 
the basic (complex) multipoles C a , E a and M a : 

Wj ( Pl = E^W $ ^(o o)&-i>.J-i I &* ( 4 °) 
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W?£ = T, P &' + J'*°*( J U' v,j (41) 



W TL{+) 
VV JJ'L 



= -T, P i + + i' + J^'A i -i o l^-r.J-^ + G-r.j-^) (43) 



(7 CT 



w£fa> = wjw (44) 

= (45) 



a' a ^ ' 

x (^-r^-iC - G-*',7-^v) (46) 



x(G-^-^v-G-^-^Jv) (47) 
^ - ^ (48) 
= Wj^ (49) 



-tt(-) _ _v-p+ , ( 3 J ' L 

JJ'L — 7 J r M'+J'\ 1 > I l l -2 

a' a ^ 

X i^J'-l'J-l^a'a ~ ^J'-l'J-l^cr'a )■ (^1) 



2.3 Angular reduced response functions 

The (q, w)-dependent reduced response functions discussed above completely determine the 
responses; however, it is not always very convenient to work with them because their number 
is infinite. In fact, the multipole expansion in Eq. (6) of the nuclear wave function involves 
an infinite number of partial waves I. There are interferences between different partial 
waves I and /' in the response functions, where I and V are effectively coupled to total 
angular momentum J' . As a consequence, the sum over J' in Eqs. (34-39) is also infinite. 
The same happens for the photon angular momenta, labeled by the multipole values J, J' 
coupled to L, that are also summed up to infinity. In practice, of course, we actually sum a 
(usually) large but finite number of terms until convergence is reached. 

For experimental separation purposes and also for theoretical analysis it is convenient to 
define an alternative finite set of what we call angular reduced response functions using the 
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fact that the values of J run from J = 0, . . . , 2 Ji. Then an expansion in Legendre functions 
Pj 4 (cos 9* ) and trigonometric functions of the angles (f>* and <p' has a finite well-defined 
number of terms, from which a useful set of response functions (coefficients of the expansion) 
that depend on q, w and the emission angle 9' can be defined as some combinations of the 
functions Wj^l(q,uj) and Legendre functions of the angle 9'. To do this we note that all 
of the angular dependence of the response functions is through the coupling of two spherical 
harmonics that can be written as 

[Yj(n*)Yj,(p')] LM = £ (JMJ'M'\LM}Y JM (e*,0)Yj, M ,(e',0)e iM *'- iMA 't', (52) 

MM' 

where A(f> = (f>' — (f>* is the azimuthal emission angle measured with respect to the polarization 
projection on the scattering plane. Note that the spherical harmonic for azimuthal angle 
equal to zero, Yjm (9*, 0), is proportional to a Legendre function Pj 4 ^(cos9*). If we take 
the real and imaginary parts of the above functions, as defined in Eq. (32), we have 

3 

Ajj, LM (n*,p') = cosMcj)' J2 h J J' lm(0')Pj (cos 9*) cos(MA^ 

M=0 
3 

+ sinM<// h^ i j, LM (9')P^ 4 (cos9*)sin(MA(j)) (53) 

M=0 
3 

Bjj, LM (n*,p') = sinM0' h^j, LM (9')P^ (cos 9*) cosiMAcj)) 

M=a 
3 

- cos M<j)' hJ 33'LM (O')Pj 1 ( cos 9* )sm(MA(t>). (54) 
M=a 

These equations can also be taken as the definitions of the angular functions hj t j, LM (9') 

and h J jj, LM (9'). Inserting Eqs. (53,54) into Eqs. (34-39) it is straightforward to write the 
response functions in the following way: 

TZ L = W L (55) 

K T = W T (56) 

K TL = cos (t)'W TL +sm<j>'W TL (57) 

TZ TT = cos 2(f>' W TT + sm2(f>'W TT (58) 

U T ' = W T ' (59) 

TZ TL ' = sm<P'W TL ' + cos <p'W TL ' . (60) 

The nine response functions W K and W K without any sub-index depend on q, ui, 9' , 9*, 
and Acj). Their explicit dependences on the angles (9* , A</>) are the following: 

2J, 3 

WL = ^T,T,fj P J 4 ( cosd ^[ P 3 c MW^ ) +P jSM W^] (61) 

J=0M=0 
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2J t J 



W T 


= ^ E fjPPicose*) [p+cmW^ + Pj SM W T jti] (62) 

J=OM=0 


W TL 


= 47r£/>P^(cos0*) 

JM 


[PjCmW^ + PjSmW^- 


)" 


(63) 


W TL 


JM 


'PpMW^ +) -PjC M W T j%- 


) 


(64) 


W TT 


= 47T^/>P^(cO S r) 
JM 


'p+cmW^+PjSmW^ 


) 


(65) 


W TT 


JM 


[p+smW^-P-cmW^- 


) 


(66) 


W T ' 


= 47T^/>P^(cO S r) 
JM 


'PjC M W^+P+s M W^[ 




(67) 


W TL ' 


= ^y^PjPj 1 ^*) 

JM 


J 3 M vv JM r J C M vv JM 


' (68) 


W TL ' 


= 47T^/>P^(cO S r) 


p-r^W TL '^ ) 4- P+^,W TL ' {+) 

t-j c M yy j M -r tjSm yy jm 


(69) 



JM 



where we define sm = sm(MAcj)) and cm = cos(.MA0), the factors that contain all of 
the dependence on A<j>, and where for brevity we have only written the ranges of the sums 
over J, M. in the L and T pieces. In all of the others a sum tov^J = 0, . . . , 2 Ji and for 
non-negative values of M. must be understood. The W's and W's inside the sums are 
the angular reduced response functions. Note that for the limit values 6* = CP, 180° the 
Legendre function P^ /, (cos0*) is zero for M ^ 0. Therefore the functions W K , W K do not 
depend on A<j> for these values of 9* . Note also that the angular reduced response functions 
Wjj^\q,u>,0') are linear combinations of the reduced responses Wjj^ l with the functions 
h'tjj'LM as coefficients, while the angular reduced responses with a tilde Wjj^l have as 
coefficients the functions ^^j-/ iM - Furthermore, in Eqs. (61-69) the responses without tildes 
are multiplied by cos(MAcf)), while the responses with tildes are multiplied by sin(.MA0). 
They can be written as 

W^\q,w,e') = Y, a ^^ h JJ'LM(d')W^l M (q^) (70) 

J'L 

W^\ q ,u,e') = Y, a ^'^J i J-LM{0')W^P LM { q ,u), (71) 

J'L 

where the index M depends on the particular response as 

( iorK = L,T,T' 
M ~ < 1 for K = TL, TV (72) 
[ 2 for if = TT, 
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that is, the value of M depends on the particular response in the same way as in Eqs. 
(32,34-39). Finally, we have introduced a factor clj>l = Pj>+l f° r K = T, cij>l = Pj> + l 
for K = T', while aj>L = 1 for the other responses. 

From these general expressions, one can easily specialize to some important cases, such 

as: 

1. Nucleon knockout from spin-zero nuclei. In this case in the above equations we have 
Ji = 0, and thus only the terms J = 0, J' = L, Jf = J and Jf = J' contribute. For 
brevity we do not give the explicit expressions for this case, as in this work we are 
considering only odd nuclei with half-integer spin. 

2. Inclusive responses. Integrating over the emission angles p' and summing over the 
final nuclear states \B) (including over proton and neutron emission) we obtain the 
inclusive responses. Due to the presence of the spherical harmonic Yj>m'(p'), the 
integral over p' in all of the above equations restricts the values to J' = 0, J = L, 
j = j', Jf = J j, A4 = M and 1 = 1', yielding an alternative derivation of the equations 
given in Ref. [6] and also providing a test of the multipole expansion. 

3. Unpolarized responses for any spin Ji nucleus. The unpolarized responses can be 
computed as the polarized responses averaged over all polarization directions, that is 



and therefore only the terms with J = in the above equations survive. These terms 
are obviously independent of fi* as is required for an unpolarized observable, so the 
integral trivially yields a factor 4tt that cancels with the denominator in Eq. (73). As 
a consequence, the unpolarized responses are given by the terms with J = A4 = in 
Eqs. (61-69). 

3 Formalism for one-hole shell model nuclei 

As a particular application of the formalism introduced in the last section to a simple model 
of a polarized target, in this work we consider the case of a one-hole shell model nucleus. 
Specifically we consider the same model introduced in Ref. [6] for the analysis of inclusive 
polarized responses and refer the reader to the discussions presented there for more details 
on the model. 

The initial nuclear state is a hole in a closed-shell core \C): 



where \i) — \rii, \, h,ji) is a single-particle state occupied in the core and b\ m . is the creation 
operator for a hole. Here we follow the convention of using lower case letters ji for half- 
integer angular momenta. As in the present work we only consider the one-body piece of 




(73) 




(74) 
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the electromagnetic nuclear current, the interaction with the virtual photon gives rise only 
to particle-hole (p-h) excitations. Thus the final nuclear states are described by 

\f) = \(E'lj)(h-\r 1 )J B ;j f ). (75) 

Here \h) — \nh,lh,jh) is another (bound) single-particle state in the core, while \E'lj) is a 
single particle in the continuum. The residual nucleus is represented in this simple model 
as a two-hole nucleus \B) = i^ 1 ) Jb) with total angular momentum Jb, coupled with 

the outgoing particle \E'lj) to a total angular momentum jf. 

In the present work the wave functions of the single-particle hole states are obtained using 
a mean-field potential of Woods-Saxon type. More details on this aspect of the calculation, 
including the values of the potential parameters, can be found in Refs. [9] and [11]. On 
the other hand, the wave function for the ejected nucleon is obtained as the solution of the 
Schrodingcr equation for positive energy e' = E' — M {E' contains the rest mass), using a 
complex optical potential in order to include absorption processes in the FSI that are not 
described by a real mean-field potential. Below we give more details about this particular 
aspect of the modeling. 

The basic quantities required are the reduced response functions in Eqs. (40-51). The 
sums in those equations run over a = J,jf) and a' — (l',f, J',j'f), but as in the inclu- 
sive case [6], we are able to perform the sum over jf,j'f analytically, simplifying the final 
expressions and permitting us to draw several conclusions about the responses. We sum- 
marize the procedure to be followed in Appendix B. Worth mentioning here is the fact that 
in Ref. [6] where we considered only the inclusive case we were able to perform analytically 
an additional sum over the daughter angular momentum Jb, considerably simplifying the 
procedures: here we cannot perform this additional summation since the value of J B is fixed. 

It is noteworthy that in this model of coincidence scattering, after performing the sums 
that are permitted, the resulting nuclear responses can be expressed in terms of the polarized 
response functions for a single particle in the shell h. These single-particle responses are 
due to the scattering of the electron with a single nucleon in the shell h, and in this case 
Eq. (19) is also valid, with the exception that Ji = jh and J B = 0, as explained in Appendix 
B. We denote the single-particle reponses of the shell h by u>j^[h — > 0]. First we consider 
the case h ^ i. Then, as shown in Appendix B, the corresponding response for the entire 
nucleus, that we denote by Wjffi [i^ 1 — > (h^ 1 ^ 1 ) Jb] , is proportional to the response of 
a single nucleon in the shell h via the simple relationship 

wfff [r 1 - (h-H- 1 )^] = (-i)^+^[J B ] 2 { I I j Jb } w^[h _ 0]. (76) 

Actually in Appendix B we prove the result at the more elementary level, namely between 
the reduced responses Wjj^l [i^ 1 — > (/i _1 j _1 ) Jb] and Wjjj\[h — > 0]. 

This relationship has important consequences. Indeed, although in this model the po- 
larization of the nucleus is carried by the outer shell i, it is interesting to note that nucleon 
knockout can provide information about the polarized responses of the inner, complete (un- 
polarized) shell h, as long as the 6-j coefficient is nonzero. In particular, if i is a neutron (or 
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a proton) we can measure polarization observables with (e, e'p) (or (e,e'n)) reactions. The 
reason for that behaviour is that the final particles are partially polarized because they are 
coupled to a definite angular momentum Jb- 

For further insight into this last comment, let us consider the sum over all possible values 
of Jb in Eq. (76). Then we obtain the sum rule for the angular reduced responses 

E W JM [I' 1 - {h-H^)J B } = 5 J0 5 M olMj h }w« (±) lh -+ 0], (77) 

Jb 

and a similar result for the total responses. For instance, using Eq. (61) for the longitudinal 
case we have 

5> L [r 1 (h-^Js] = 4*\j h ]w$ +) [h 0], (78) 

Jb 

where we have used /q = I /[.]%}■ We can compare this equation with the longitudinal 
response for electron scattering from an unpolarized particle in shell h. From Eq. (73) and 
the comments that follow it, we see that the unpolarized response is given by the term with 
J = 0in Eq. (61), with Ji =j h 

Knpol [h^0}= 4^«4 (+) [h^O]. (79) 

Then, inserting this equation in Eq. (78), we obtain a relation that is valid for all of the 
responses, K = L,T, . . ., and not only for the longitudinal case: 

J2k k [r 1 - [h-h-^Js] = lj h } 2 n« ipol {h -> o]. (so) 

Jb 

This shows that the sum over Jb of all of the exclusive polarized i^-responscs for different 
final states of the type (h~ 1 i~ 1 )JB is equal to the unpolarized ^-response of the shell h 
(that is, arising from unpolarized nucleon knockout from a nucleon in the shell h) times the 
number of particles in the complete shell h. As a consequence, the sum over Jb of all of 
the polarized responses (i.e. with J ^ 0) is equal to zero. Of course this theorem is only 
valid in the extreme shell model, which means that observable deviations of this result can 
provide valuable information about configuration mixing and the nuclear structure of the 
final states. In summary, we see that the reason why unpolarized shells yield polarization 
observables is due to the fact that the angular momentum Jb of the daughter nucleus is 
fixed. If we sum over Jb all of the polarization observables go away for a closed shell. 

Now we consider the case h = i, where we extract a nucleon from the same shell in which 
the hole is located. This is the case with more practical applications, because the shell model 
is expected to provide a better description of the nuclear states for small excitation energy 
and, in particular, for the ground state of the daughter nucleus with Jb = 0. The derivation 
of the results that follows is similar to the h ^ i case, with the exception that a factor of 
two arises due to the fact that the two holes in the same shell are undistinguishable, and, in 
addition, that the nuclear spin Jb is restricted to be an even number. Then, if we set h = i 
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in Eq. (76) and multiply by two, we obtain 

W^[r^(rH-i)J B ]=2(-l)^[J B f{ * j f U^[^0]. (81) 

In particular, for the daughter nucleus in the ground state with J B = the 6-j symbol is 
equal to (— l) J+1 /[ji] 2 and we have 

Wfff [r 1 - (r^jO] = ^^[i - 0]. (82) 

For example, for the simplest case of a hole in a ji — \ shell, we have [ji] 2 = 2 and therefore 
we obtain 



W K{±) 



2 



(rV)° =^S } [i-o]- (83) 



Then we recover the result that one hole in a i shell is equivalent to a particle in that shell, 
as expected. 

Also, if we sum over J B as before, we obtain an interesting sum rule. Taking into account 
that the sum now runs over J B = even and using the property 

j B =even k ' 

then, from Eq. (81), we obtain for the sum of the angular reduced responses 

E W JM - ii-^Js] = (6jo\ji? (-l) J ) W f^[ t - 0], (85) 
Jb 

i.e., the unpolarized response of the 2j t + l particles in the complete shell minus the polarized 
response of a single particle multiplied by (— 1) J . The same relation holds between the sum 
^2 j Wjjf^ [i^ 1 — > (i~ 1 i~ 1 )J B ] and the responses WjjJ^i — > 0]. The origin of the phase 
(— 1) J is easily explained by looking at the total responses. For instance, from Eq. (34) we 
obtain for the longitudinal one 

[r 1 - (rt- 1 )./*] = \ji] 2 K npol [i^o] 



Jb 



-47T E fj(-V J {PjAjJ'MVjpLli ^ 0] - PjBjj, LO w L j { -) L [i -+ 0]} . (86) 

JJ'L 

We notice that we can include the factor (— 1) J in the angular functions Ajj'lo and Bjjilq 
using the parity property of the spherical harmonics Yjj^(—ft*) = (— 1) j Yjm(Q*)- Then 
we can write the above result for any response, and not only for the longitudinal 

[r 1 - [i-H-')j B ] = \ji] 2 nZ npol [i - o] - - o]| n ._ n . . (87) 

Jb 
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Now it becomes evident that the factor (— 1) J+1 for the difference between particle and hole 
cases which was also mentioned in Rcf. [6] comes from the fact that a hole nucleus which is 
polarized in the direction CI* is physically equivalent to the absence of a particle which is 
polarized in the opposite direction — CI* . The factor (— 1) ^ comes from a spherical harmonic 
Yj_m(Q*) which is always in front of the responses of rank J in the multipole expansion. 

A similar result to Eq. (87) was obtained also in the inclusive case (see Ref. [6]), where 
it was shown that any polarized (inclusive) response of a hole nucleus can be written as 
the unpolarized response of the complete shell minus the polarized response of the hole 
(multiplied again by the factor We can re-obtain here that result by integration of 

Eq. (87), thereby showing another connection between the two formalisms. 

To finish this section we show a result for the sum of the unpolarized (exclusive) responses 
of a hole nucleus: 

J2 U un P ol - (r'rVfl] = m 2 Wun P oS - 0], (88) 

Jb 

i.e., the sum of the unpolarized responses for each one of the particles in the shell. 

Equations (76) and (81) are the basic results of this section. They show that the reduced 
polarized responses corresponding to different values of the final angular momenta of the 
daughter nucleus \(h~ 1 i~ 1 )Js) are scaled by a well-defined coupling coefficient that only 
depends on the values of the angular momenta involved. This fact could be experimentally 
checked and any deviation from that relationship would imply interesting conclusions about 
the nuclear structure of the daughter states. 

4 PWIA 

4.1 Factorization of the responses 

The main reason for considering the PWIA for the ejectile is that in such a limit, when the 
final-state interaction is turned off, the exclusive cross section factorizes and we can perform 
part of the calculations analytically. The resulting cross section contains the basic electron- 
nucleon (in general off-shell) responses and the scalar and vector momentum distributions, 
giving a picture of the physics that is more transparent than the general formalism where 
the multipole analysis obscures the physical interpretation of the process. 

In order to extract valuable information about the nuclear momentum distribution, it 
would be desirable to find a set of polarization observables that is more or less independent 
of the FSI. The general formalism for PWIA in polarized nuclei was developed in Refs. [5, 7], 
and it was recently specialized in Ref. [6] to the particular case of polarized, one-hole nuclei, 
although in that work the exclusive responses were integrated over the emission angle in 
order to study the inclusive case. Here we briefly summarize the essential results from that 
work so that the nomenclature can be introduced, referring to Ref. [6] for more detailed 
discussions. 

In PWIA we write the factorized exclusive responses as 

TZ K = Mp'w§M s (p) (89) 
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K K ' = Mp'wf -M y (p), (90) 

where p = p' — q is the missing momentum with angles (6, <fi), and Wg and Wy are respec- 
tively the scalar and vector single-nucleon responses, which are related with the components 
Eqs. (9-14), of the single-nucleon tensor spin matrix 

^(p',p) = ^(p',p) t ^(p',p), (91) 

where J M (p',p) is the electromagnetic current spin matrix for the nucleon, evaluated with 
plane waves (r|ps) = (27r)~ 3//2 e 4p ' r x s - Taking the appropriate components of the single- 
nucleon tensor w' il/ (p',p), the resulting unprimed responses are diagonal in spin (spin- 
scalar), so we can write for the L, T, TL and TT responses w K = Wg , while the T 1 and 
TV (primed) responses are spin- vectors and can be written as w K — Wy ■ cr. Explicit 
expressions for the current matrices and single-nucleon responses are given in Ref. [6] . 

The momentum distribution for residual nucleus \B) is defined in general as the spin 
matrix 

n(pVa - Y,(B\a p , s \A)*(B\a p , s ,\A), (92) 

M B 

where the spin Jb is fixed and we sum over undetected orientations of the daughter nucleus. 
Therefore the definition of the scalar M and vector momentum distributions is 

n(p) = ^(M s (p)+M v (p)-a). (93) 

Obviously, the polarized momentum distribution depends upon the initial state \A), the final 
state \B) and the polarization angles ft*, apart from the missing momentum p. When wc 
need to clarify this dependence below, we sometimes write in the momentum distribution 
the explicit process to which we are referring as an argument, for example, n [A(Q*) — ► B]. 

4.2 Momentum distributions of hole nuclei 

Here we obtain the polarized momentum distribution in the particular case of a hole nu- 
cleus. For the general expression of the polarized spectral function of deformed nuclei, 
see Ref. [7]. In this work we are also interested in examining in more detail the relations 
between the single-particle momentum distribution and that of the whole nucleus in the 
extreme shell model, and thus we will follow a somewhat different approach (although one 
that is completely equivalent). For this reason we begin examining the polarized momentum 
distribution of a single particle (or hole) in a shell. 

4.2.1 Single-particle momentum distribution 

Let us consider a single particle in the bound state = J2 m ^mjb m ) °f a mean field; 

here one should understand there to be the argument fT in all the rotation matrices, polar- 
ized states and momentum distributions. In this simple case, the momentum distribution 
can be written as the matrix element 

n s , s \j(Q*)} = (ps'\N\j(n*)]\ps) (94) 
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of the spin density operator N[j(Sl*)] which projects onto the state \jj*) 

N\j(Q*)\ = \m)(n * I = E ^^S-WXiH- (95) 

mm' 

Expanding the product of the two rotation matrices, we are able to write a multipolc expan- 
sion for the operator N[j(il*)] in a basis of spherical harmonics of the polarization angle: 

N\j(n*)} = J2fj Nj Un (96) 
J 

where N J {j(Q*)] is the operator of rank J in the expansion, given by 

N*m*)\ = yfo E (-l) 3+m+J ( _! m 3 m , J M ) Y JM (Q*)\jm')(jm\. (97) 

mm' M 

A similar expansion can be written for the case of a single hole, namely, the focus in 
the present work. In the shell model a hole is described by a single-particle state with wave 
function \jj*) = J2 m ^rnj f \j m )i i- e -> the time inversion of the single-particle state, with 
\jm) — (— iy +m \j, — m). In this case one can prove that the corresponding spin-density 
matrix operator multipole expansion can be written as 

N\j(fc)] =J2fj(-l) J N*m*)}, (98) 
J 

were we have used the fact that the spin-density multipolcs of particles and holes differ only 
by the phase factor (— 1) J . This a consequence of the fact that the state |j, —j) rotated by 
an angle SI* is the same as the state \jj) rotated by an angle — Q* and of the parity (— 1) J 
of the spherical harmonics which are the basis of the multipolc expansion. 

4.2.2 Hole nucleus 

In analogy with the single-particle case, we can write the momentum distribution as the 
matrix element of a one-body density matrix. In fact, consider the transition 

\A) = r\W)) ^\B) = \{h-\^)J B ). (99) 

Inserting the appropriate matrix elements of the destruction operator into the momentum 
distribution, Eq. (92), we obtain 

n s , s ( P ) - ( P *'|jv - (/^r 1 )^] |p S ), (loo) 

where N [i _1 (0*) — > (h~ 1 i~ 1 )Js] is the spin density matrix operator 

wfr 1 ^*)-^- 1 *- 1 )^] = 

E E E (3hm h j i m^J B M B )(j h m' h j i m^\J B M B )V^?*V^ (101) 

Mb wthm| m' h m'f 
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Expanding again the product of the two rotation matrices we obtain a sum of the product 
of three 3-j coefficients, that can be written as a 6-j coefficient times a 3-j coefficient. As 
consequence, we straightforwardly obtain the multipole expansion 

n [r\n*) - (h-'t- 1 ) j b ] = Y, fjN J - {h-H- 1 ) j b ] , (102) 

j 

with the simple relationship between the nuclear multipoles and the single-hole h results: 
N J -» (h-H-^Js] = [J B ] 2 (-l) ji+ih+JB+J \ jh jh f }N J [h(Cl*)}. (103) 

This relation is the same as the one found between the reduced nuclear responses and the 
responses of the shell h (see Eq. (76)). In fact, the general equations of the last section are 
of course also applicable to the PW1A when we take the single-particle state \E'lj) to be a 
free wave function, i.e, proportional to a Bessel function. Therefore, Eq. (76) had to be valid 
also in PWIA. And that fact is again recovered here in another completely different way, 
providing a test of both formalisms. But also the fact that the relation between momentum 
distributions can be extended to response functions even with FSI leads us to believe that 
some information about the single-particle momentum distribution can be extracted from 
these kinds of observables. It is remarkable that, although the FSI involving the ejected 
nucleon destroys the factorization of the cross section, it does not modify the relation in 
Eq. (103), implying that this property is general for the independent-particle model. 

Until now we have only considered the case h =/= i. The other situation h = i that arises 
when the particle is ejected from the polarized shell is analogous, with the exception that 
now the residual nucleus has a wave function with a factor 

\B) = -j=[blbi]j B M B \C). (104) 

Therefore the expression for the momentum distribution as the matrix element of a spin 
density matrix operator is 

n s , s [r x (o*) -» (i-H-^b] = ( P s'\N [r 1 ^*) -» [rH-^B] \ ps ) (105) 

with 

N J - (i-H- 1 ) J B ] = 2[J B ] 2 (-1) 1+J ( k h f } N J [i(n*)}, (106) 

and we again recover the relation in Eq. (81) between reduced response functions in the shell 
model. 

These relations between the multipole components of the density matrix lead us to obtain 
sum rules when summing over the the values of Jb , as we did for the responses in the shell 
model. In fact, performing the sum over Jb in the case h ^ i we obtain a the result analogous 
to Eq. (80) 

-> (h-H-^Jg] = [jh} 2 N[h] unpol , (107) 

Jb 



19 



where N[h] unp oi in the umpolarizcd density matrix of the single-particle wave function. This 
result was also obtained in Ref. [6] in relation with the inclusive PWIA, where the final state 
B is not observed and the sum over Jb must be performed. Strictly speaking, what in Ref. [6] 
is called polarized momentum distribution 

of a shell actually is the matrix element of the sum over Jb given by Eq. (107). Then 
we can say that the momentum distribution of a complete shell (h ^ i) is the sum over all 
of the unpolarized momentum distributions of all of the particles in the shell. 

In the case h = i we have the restriction J B = even. The sum over J B gives an equation 
similar to Eq. (85), and the total density matrix is 

Y, N [r 1 ^*) [h-H-^JB] = \ji] 2 N[i\ unpol - N[i(-n*)), (108) 

Jb— even 

which corresponds to Eq. (87) at the response level. This result was also obtained in Ref. [6] 
for the polarized momentum distribution of a shell with a hole: it is equal to the unpolarized 
momentum distribution of the complete shell minus the polarized momentum distribution of 
the hole, although in this work we are able to explain geometrically the origin of the factors 
(— l) J that appeared in that reference and also to understand the fact that the momentum 
distribution of a hole is equal to the momentum distribution of a particle pointing in the 
opposite direction. 

4.3 Scalar and vector momentum distributions 

Until now we have seen that when the polarized momentum distribution (or in general the 
spin density matrix operator) of hole nuclei is expanded in multipoles of the polarization 
angles, each one of the individual multipoles is proportional to the corresponding multipolc 
of the single-particle momentum distribution of the struck nucleon. Therefore the same 
proportionality relation also holds for the scalar and vector pieces in Eq. (93) of that spin 
operator. Thus we only need to describe the scalar and vector momentum distributions 
of a single particle in the state A detailed derivation of these expressions was 

performed in Ref. [6]. Therefore here we only write the results in a way more suited to use 
in the exclusive responses. The multipole expansion for the scalar and vector momentum 
distributions reads 

M s [j(Q*)} = Yfj P J M jm*)\ (109) 
J 

J 

where in addition we have used the result that only even (odd) multipoles enter in the 
expansion of the scalar (vector) momentum distribution for a single particle. Although not 
explicitly written, the functions Mj and M \ depend on (fi* , p) . The only scalar that can 
be constructed with Yjj^(Q*) and p is [Yj(Q*)Yj (p)]oo, therefore the scalar momentum 
distribution for a single particle has the form 

M|(fi*,p) = m S j(p)[Yj(n*)Yj(p)} 00 . (in) 



20 



Also, the vector momentum distribution can be written as a linear combination 

M^(Q*,p) = rn^j,(p)Xjj,(n*,p) (112) 

J'=j±i 

of the vectors ~Kjj'(£l* ,p) with spherical components 

Xjj,(<n*, P ) a = [Yj(n*)Yj,(p)] la . (113) 

Finally, the functions m.j (p) and nij j, (p) for a particle are given by 

m S j\j{W)\ = 2Aj[j] 2 \R(p)\ 2 (H4) 
m V jj,m*)] = -^Ajj,[jf\R(p)\^ (115) 

where R(p) = \f2~pn dr r 2 ji (pr)R(r) is the radial wave function in momentum space, 
normalized to dpp 2 \R(p)\ 2 = 1, and the factors Aj and Ajj> are coupling coefficients 
which can be found in Ref. [6]. Finally, the same relations in Eqs. (111,112) are valid for 
the scalar and vector momentum distributions for the process i _1 (J7*) — > (h~ 1 i~ 1 )Js and 
we have the following result: 

{l + 5 hi )[JB?{-V ii+jh+JB+J { i h j - f )m s j[h(n*)] (116) 

and a similar result for rrijj, [i _1 (f2*) — > (/i _1 i _1 ) Js] in terms of the single particle vector 
multipole rrijj, [h($l*)]. 

4.4 Single-nucleon current and exclusive nuclear responses in PWIA 

For the single-nucleon current we follow the formalism of Refs. [9, 6], where the on-shell 
rclativistic electromagnetic single-nucleon current was expanded up to order rj = p/M , but 
keeping all orders and not expanding in the dimensionless variables k = q/2M and A = 
lo/2M. The expressions for the charge and transverse current components in the coordinate 
system used in this work are given in Ref. [6], where also can be found the expressions for 
the single-nucleon responses w K . The important point here, and used for the results of next 
section, is that the L, T, TL and TT single-nucleon responses are spin- independent, while 
the T' and TL' ones involve linear combinations of the Pauli spin matrices. 

We now have all of the ingredients needed to obtain the reduced responses of a one- 
hole nucleus in PWIA. The following equations have been derived for a one-hole nucleus, 
although actually they are a little more general — in fact, they are valid for any nuclear 
model as long as the multipoles of the scalar and vector momentum distributions can be 
written as in Eqs. (111,112). 

For comparison with the general model it is convenient to write the responses in a form 
similar to Eqs. (61-69). Therefore, for the unprimed responses that are proportional to the 
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scalar momentum distribution M s , we develop the (real) scalar coupling of two spherical 
harmonics as in Eq. (53): 

J 

[Yj(Q*)Yj(p)} 00 = J2 ^(cosr)cos[MA0]^ oo (0), (117) 

M=0 

and so we have 

J 

Mf(fT,p) = ]T P^(cose*)cos[MA^h^ J00 (e)m S j(p), (118) 

M=0 

where rrij refers to the considered transition (see Eq. (116)). For the vector momentum 
distribution we need to develop the three components of the vector ~K.jj> in terms of the 
functions h-jj,^, h-jj,^ and hjj, w in a similar way. 

As for the nuclear responses, using Eq. (118) we see that the scalar momentum distribu- 
tion depends on the azimuthal angles only via cos MA<f), where A<f) = (f)' — <fi* , but not on 
cf> = (f>' or (f>* independently. The only dependence of the nuclear responses on <p is therefore 
through the single- nuclcon responses w TL and w TT , see Ref. [6]. Then it is straightforward 
to read off the following expressions for the angular reduced responses: 

Wffl = ^(Pl+Pls^j^m^p) (119) 

Wjtf = ^(2Jl + J^)hy jm (e)rn S j(p) (120) 

Wj L A +) = ^-2V2(p c J c + p so J m )Sh^ J0Q (e)m S j(p) (121) 



47T 

- -^J^h-^H^^. (122, 



r TT{+) _ Mp' T 2 x 2 h M^ im .„S, 

7TL{+) , wTT(+) 



The odd-rank unprimcd responses and the Wj^ ' and WjjJi ' responses are zero in PWIA. 
Here 5 = 77 sin labels the order of the relativistic correction and (9,<f>) is the direction of 
the missing momentum p. The factors p c (charge), p so (spin-orbit), J m (magnetization) and 
J c (convection) come from the electromagnetic charge and current operators; they are only 
(<7,a>)-dependent, and are defined by (see Ref. [6]) 

Pc = —j=Ge (123) 

Jm = V^Gm (125) 
Jc = —G E , (126) 

K 

where Ge, Gm are the electric and magnetic form factors of the nucleon for which we use 
the Galster parametrization [12]. 
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Finally, for the primed responses, after performing the scalar product between M v and 
Wy we obtain 

W T jtf = Y, [jlhfj.^ + ^JMhfj.^rn^Ap) (127) 

J'=J±i 

WjM { ~ ] = -^f E ^PcJ^,^ +PsoJ m 5hy j , w {6)]ml J ,(p) (12%) 
V2tt j/=j±1 L J 

= E (PcJ m -p ao Jc« 2 )^nW"»^(p). (129) 

All of the other responses not written here are exactly zero in PWIA. 

5 Results and discussion 

In this section we present numerical results for the reaction 39 K(e, e'p) 38 Ar in the shell 
model for both PWIA and with FSI, namely the distorted-wave impulse approximation 
(DWIA). An overview is provided in Sec. 5.1, followed by Sec. 5.2 where detailed results are 
shown and discussed. Section 5.3 contains a brief discussion of the relativistic corrections 
incorporated in our model. Throughout all of these sections we have had to select only 
the most important results for presentation, since the entire set vastly exceeds what can be 
shown in published form. A more complete set of figures may be found on the WWW [21]. 
The reader is directed to the conclusions in the next section for discussion of what can be 
learned from such studies. 

The bound states of the target and daughter nuclei are represented by one-hole and 
two-holes in the core of 40 Ca, as discussed in Sec. 2 and in Ref. [6]. It is clear that more 
sophisticated treatments of nuclei such as 39 K and 38 Ar than the extreme shell model we are 
using in this work can also be studied, however, since our aims here and in the next paper [19] 
are to analyze how FSI effects are reflected in the various polarization observables and to 
determine if it is possible to obtain information about the different pieces of the polarized 
momentum distribution or about the FSI with these kinds of measurements, this basic 
model should provide an adequate starting point. Indeed, given the complexity that already 
exists for our basic model and the fact that no such extensive study for a complex nucleus 
is available, the initial studies presented in Refs. [7, 15, 8] being all that we are aware of 
for coincidence polarization observables, at this time explorations with more sophisticated 
nuclear configurations are probably unwarranted. Of course, when experimental studies of 
this type become imminent these extensions can be undertaken. 

In our modeling the single-particle bound wave functions and (negative) energies are 
obtained using a mean-field potential of Woods-Saxon type, using the values of the potential 
parameters that may be found in Refs. [11, 16]. For the ejectile wave function we solve the 
Schrodinger equation for positive energies using a (local) complex optical potential fitted 
to proton elastic scattering data. For the particular case of the nucleus 39 K we use the 
potential of Schwandt et al. [17], that was fitted to cross section and analyzing power 
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angular distributions for a number of target nuclei and polarized protons of 80-180 MeV. 
Although the daughter nucleus 38 Ar was not included in this analysis, it falls in the range 
24 < A < 208 of applicability of the Schwandt potential. As a first test of our model, we 
have computed the cross section and analyzing powers for proton scattering from a variety 
of nuclei with different optical potentials and have compared the results with the literature. 
In particular, we have reproduced the calculations shown in Figs. 2 and 3 of Ref. [17] for 
the analyzing power of protons clastically scattered from 24 Mg, 28 Si, 40 Ca, 90 Zr and 208 Pb 
at different energies, using the Schwandt potential. 

The normalization of the partial wave with energy e and wave number k — \/2Me 
is determined by the following asymptotic condition for the radial wave function 

R lj( r ) ~ \ ^-c^'+ s '^ sin (kr -f]\o g 2kr - I J + a t + , (130) 
V irh k ^ 2 / 

where Sij is the complex nuclear phase-shift and m is the Coulomb phase-shift. In the 
limit when the imaginary part of the potential goes to zero, the radial wave functions are 
normalized with a Dirac <5-function containing the energies, so that condition in Eq. (8) is 
satisfied. The imaginary part of the potential modifies the normalization of the continuum 
states, since the flux of outgoing particles is reduced due to the absorptive property of the 
optical potential. In addition, in this model the continuum nuclear wave functions are no 
longer orthogonal to the bound ones. However, at high momentum transfer q <~ 500 MeV/c 
and near quasi-free conditions where p' ~ q, the orthogonality between the bound and 
high-momentum ejected particles is approximately satisfied. Several prescriptions to restore 
orthogonality and to take into account non-locality effects have been proposed, such as the 
inclusion of Perey-like factors [18] in the wave functions. In our calculations we include a 
Perey factor of the kind 

-1/2 

(131) 

where Vc is the central part of the optical potential and we take the non-locality parameter 
= 0.85 fm. 

In order to take into account relativistic effects that are important for high momentum 
transfer we use relativistic kinematics in all of the calculations [6]. For nuclcon knockout 
from shell h we compute the momentum of the ejected particle with kinetic energy e' = eh+uj 
using the relativistic energy-momentum relation p' 2 = (M + e') 2 — M 2 . In the shell model 
we solve the Schrodingcr equation with eigenvalue e'(l + e'/2M). 

In DWIA the sums over partial waves or multipoles J are infinite, while in practice we 
must sum all multipoles until convergence is reached. In this work we test the degree of 
convergence by comparing with the exact result in factorized PWIA. In fact, in the limit 
in which the optical potential goes to zero, the DWIA and PWIA approaches should be 
identical and therefore in that limit we adjust the number of multipoles so that differences 
with exact PWIA responses are indistinguishable. The number of multipoles needed to 
reach convergence is in general greater that in the inclusive case [6]. For example, for 39 K 
at q = 500 and 700 MeV/c in the region of the quasielastic peak it is enough to include 
multipoles up to J = 15 and 23, respectively, for the inclusive responses, while for the 



f(r) 
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exclusive process we need to sum up to J = 22 and 32, respectively. The reason is that in 
the inclusive case the integration over the angles p' of the ejectile is performed analytically. 
In all of the calculations that follow we include multipoles up to J = 32. 

5.1 Overview of the responses 

We apply the formalism of previous sections to the nucleus 39 K. In our model we describe 
the ground state as a proton hole in the shell ld 3 / 2 . Therefore we have Ji = 3/2. Since 

J = 0, . . . ,2Ji, the angular reduced response functions Wjffi and Wjffi that enter in 
Eqs. (61,69) range from (J,M) = (0,0) to {J,M) = (3,3). Taking into account that all 
of the responses with a tilde are zero for M. = (because they are defined as the coefficient 
of sin.MA</>), we are confronted with a total of 72 angular reduced response functions (only 
32 of them are nonzero in PWIA). 

In this work we focus our study of polarization observables upon a systematic analysis 
of the complete set of angular reduced response functions, which contain all of the informa- 
tion about the process. An analysis of more direct observables, such as cross sections and 
asymmetries for different values of the polarization angle is in progress [19]. 

We discuss only the case in which the daughter nucleus 38 Ar is in the ground state 
Jb = + , which we describe as a two- hole nucleus. Under knockout from the 1^3/2 
shell the daughter nucleus also can be in a Jb = 2+ state, but as proven in Sec. 3, the 
reduced response functions Wjff for the two possible values of Jb differ only in a constant 
recoupling coefficient (see Eq. (81)) and therefore they are proportional in this simple model, 
so no new information is obtained from their analysis. Of course, in a more sophisticated 
treatment of the nuclear structure, configuration mixing will produce deviations from this 
proportionality theorem, which may prove to be useful in obtaining new information about 
the nuclear structure involved. 

In theory the angular reduced response functions Wjff and Wjff could be extracted 
by performing a number of measurements with different electron helicity h, target polariza- 
tion ft*, azimuthal emission angle <f>, and, additionally, a Roscnbluth analysis in order to 
separate the longitudinal from the transverse responses. The number of ways to carry out 
the extraction is not unique. For example, a discussion of the minimal asymmetry measure- 
ments needed to extract the reduced response functions in deuteron electrodisintegration is 
given in Ref. [20], although in the 39 K case the presence of the octupole terms J — 3, which 
do not appear in the deuteron case, modifies the extraction rules given in that reference. 
Rather than discussing the separation procedure or even to attempt to discuss the entire set 
of reduced responses in great detail, our goal in the following analysis in the present work 
is to take into account that some of these reduced responses are very small relative to the 
unpolarized ones, and hence to focus on those with the largest magnitude. We postpone 
any detailed discussion of cross sections and asymmetries for future presentation [19]. 

The angular reduced response functions are defined in Eqs. (61-69). Those without a 
tilde Wj M are multiplied in the cross section by cosMAcj), while the responses with a tilde 
Wj M are multiplied by sinA^ A</>. Therefore the latter are zero for A4 = or A<fr = 0, 
and hence the tilde responses do not contribute to the cross section when the target is 
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Ia 

q = 300 MeV/c 
lu = 40 MeV 
p' = 244.3 MeV/c 
e' = 31.3 MeV 


lb 

q = 300 McV/c 
w = 56 MeV 
p' = 301.6 MeV/c 
e' = 47.3 MeV 


Ic 

q = 300 McV/c 
cj = 80 MeV 
p' = 372.6 MeV/c 
e' = 71.3 MeV 


Ha 

q = 500 MeV/c 
cj = 110 MeV 
p' = 447.6 MeV/c 
(' = 101.3 MeV 


lib 

q = 500 MeV/c 
w = 133.5 MeV 
p' = 499.7 MeV/c 
e' = 124.8 MeV 


lie 

q = 500 MeV/c 
w = 160 MeV 
p' = 553.9 MeV/c 
e' = 151.3 MeV 


Ilia 

q = 700 MeV/c 

= 180 MeV 
p' = 592.3 MeV/c 
e' = 171.3 MeV 


111b 

9 = 700 McV/c 
w = 241 MeV 
p' = 699.9 MeV/c 
e' = 232.3 MeV 


IIIc 

q = 700 MeV/c 
u = 300 MeV 
p' = 794.6 MeV/c 
e' = 291.2 MeV 



Table 1: Values of the momentum and energy transfer (q,ui), and of the momentum p' and 
kinetic energy e' for which the reduced response functions have been computed. 

longitudinally polarized, i.e. in the direcction of q (L-polarized) or q, because in those 
cases cos 8* = ±1 and the Lcgcndre function Pj i (±l) = for M. ^ 0. Also they do 
not contribute when the azimuthal angles of the polarization direction <fi* and the emission 
direction <f> are the same (A</> = 0). As noted above, the total number of angular reduced 
response functions for nucleon knockout from the d 3 / 2 shell of 39 K is 72, which is too large 
for a complete presentation and therefore in the following we only show figures for the most 
relevant. The complete set of figures is presently available from the authors [21]. 

The polarized responses Wjffi and Wjffi depend on the three variables {q, u,p), since 
the missing energy is determined. In order to study their dependence on all of the variables, 
we have computed the entire set of responses as functions of the missing momentum p for 
nine values of (q,u>). In each one of the following figures we show a particular response 
function and each figure is composed of nine panels corresponding to the values of (q, ui) 
given in Table 1. In each one of the rows in the table the value of q is fixed to one of the 
three values q — 300 MeV/c (regions Ia, lb, Ic), q = 500 McV/c (regions Ha, lib, lie) or 
q = 700 MeV/c (regions Ilia, Illb, IIIc). We show three values of uj for each q. The second 
column corresponds approximately to the quasielastic peak (QEP) defined by p' ~ q. The 
first column corresponds to slightly below the peak, p' < q, while the third column is in the 
region above the QEP, p' > q. The values of the momentum p' and kinetic energy e' of the 
ejectile for each of the panels are also shown in Table 1. 

First notice that the Schwandt optical potential we use was originally parametrized for 
elastic scattering of protons with energies between 80 and 180 MeV, and therefore it has been 
extrapolated in some of the kinematic regions. Also the daughter nucleus 38 Ar was not used 
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in the analysis of Schwandt et al, but it is in the mass range of the nuclei originally used. 
One of the aims of this work is to analyze the sensitivity of the different responses to the FSI 
model, and, in particular, in each one of the panels in Figs. 1-15 we show four curves, each 
one corresponding to a particular model for the final proton state. Specifically, the solid 
lines have been computed with the complete optical potential, while for the dashed lines we 
have set the spin-orbit potential (SOP) equal to zero (both real and imaginary parts). Thus 
comparison between solid and dashed lines shows the effect of the SOP in the responses. 
These two curves do not include a Perey factor in the wave function, and therefore we also 
show with dot-dashed lines the same responses computed with the total potential and with 
Perey parameter (3 = 0.85 in the wave function. Finally, the dotted curves correspond to 
the PWIA, that is, they do not include FSI. 

In order to get a feeling about the relative order of magnitude of the 72 responses 
it is convenient to choose some quantity which characterizes the importance of each re- 
sponse. The total variation of a response defined as the maximum minus the minimum 
values AWjff = Wjj^\max) — W^^^min), that is always positive, is a quantity that 
is very useful for that purpose. Specifically, we take the variation of the unpolarized longi- 
tudinal AWqq + ^ as reference and therefore give the others as a percentage of this response, 
so the variation of the unpolarized longitudinal response is 100% by definition. We have 
computed the total variation of all of the responses using the complete optical potential for 
the nine kinematic conditions and the results are shown in Tables 2-4 for q = 300, 500 and 
700 MeV/c, respectively. Note that under the columns for the L, T and T" responses both 
kinds of responses (with and without tilde) are displayed because in these cases there is no 
confusion: the L and T responses without (with) tilda are always of even (odd) rank, while 
the opposite happens for the T' responses. 

With reference to the PWIA responses, only the 32 responses given in Eqs. (119-129) 
survive in that limit. Therefore, a total of 40 of the 72 responses are exclusively due to 
the FSI, and then it is expected that they will be especially sensitive to the details of the 
interaction, as we shall show below for some of the more important cases. 

5.2 Detailed discussion of selected responses 
5.2.1 L-responses 

Unpolarized response. The unpolarized response Wqq + ^ displayed in Fig. 1 serves to 
set the scale of the remaining responses. The three panels in the middle row of Table 1 at 
q = 500 MeV/c are representative and well within the range of validity of the Schwandt 
potential. The one in the middle at u — 133.5 McV (region lib) approximately corresponds 
to the QEP. The FSI (solid lines) produce a reduction of the plane-wave results (dotted 
lines) of about 35%. When we set the SOP to zero, the response increases by a few percent 
(dashed lines); also we see that the effect is bigger for higher energy. The introduction of 
a Perey factor slightly reduces the total response (dot-dashed lines), but the effect is very 
small. A significant effect of the FSI with respect to the PWIA is seen in the extreme cases, 
q = 300 MeV/c, lu = 40 MeV (region la) and q = 700 MeV/c, lu = 300 MeV(region IIIc), 
although these results should be viewed with caution, being far outside the range of validity 
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of the potential. Indeed, in the former case the real attractive part of the potential is too 
large, while in the latter the imaginary absorptive part is too large. In region Ilia (q = 700 
MeV/c, u> = 180 MeV) the spin-orbit interaction produces an increment of the response (the 
opposite to what happens in the other cases). The effect of the Perey factor decreases with 
the energy because it only depends on the real part of the potential, which also decreases 
with the energy (this conclusion is also applicable to the remaining responses). 

Vector response. The vector longitudinal response \ shown in Fig. 2, is zero in 

PWIA. Thus it is due only to the FSI and expected to be more sensitive to details of 
the potential. In addition, it is not as large as the unpolarized one, but we find that it 
is large enough to make its extraction feasible — in fact, in region lib at the maximum 
p ~ 150 MeV/c, its magnitude is about 24% of W^ + ^ (see Table 3). In regions Ib,c and 
IIIb,c, it is even bigger. For instance, we see in Tables 2 and 4 that the variation AW^~^ 
in regions lb and Illb is about 56% and 44%, respectively, of the unpolarized longitudinal. 

Without the spin-orbit interaction this response is much larger and in some cases its 
behaviour as a function of p is even different. Specifically, in the QEP region (second 
column of Fig. 2) the SOP produces a small shift to the right, which is more important for 
low q. The same happens in the region to the right of the QEP. However, more important 
effects are seen below the QEP, where the response is significantly reduced by the spin-orbit 
interaction and its shape is smoothed and broadened. The effect is more pronounced at high 
energy; for instance, at w = 180 MeV (region Ilia) the SOP produces a reduction of the 
peak of more than 60%. 

This response contributes the most to the cross section when the polarization is normal 
to q, namely at 9* = 90° because it is multiplied by a Legendre function Pi(cos6*) = sin 6**, 
and for Atfi = ±90°. Therefore we expect the spin-orbit interaction to be important under 
these conditions. For example, for coplanar kinematics (<fi = 0) and normal (N) polarized 
39 K (6>* = 0* = 90°), where A<£ = (j) - (f>* = -90°. 

Quadrupole responses. In this case the PWIA responses W^ 1 ^ are nonzero and in 
general they are large, as seen in Tables 2-4. The biggest one corresponds to M = 0, shown 
in Fig. 3, which contributes significantly for L-polarization: for instance, its size is about 
one-half of the unpolarized response in region lib (see Table 3). For low p this response is in 
general negative and here the PWIA provides insight into why this happens. From Eq. (119), 
we see that this response is proportional to the angular function /i§ 2 oo W an( ^ to the radial 
function mf (p). From its definition, eqs. (53,54) it is easy to see that the function /12200W 
is proportional to P^cosO) = i(3cos 2 # — 1), while the radial function mf (p), defined by 
Eq. (114), is proportional to the factor A 2 , given in rcf. [6], which is proportional to a 
negative 3-j coefficient. Therefore, this response includes a factor 1 — 3cos 2 6 l . Below the 
QEP we have p' < q and therefore, for p small, p is almost anti-parallel to q because 
p + q = p', and so the angle 9 is near 180° and therefore 1 — 3 cos 2 9 ~ — 2 is negative and 
the total response is negative. As p increases, the angle 9 quickly reaches a value near 90° 
degrees, and so the factor 1 — 3 cos 2 6 ~ 1 changes its sign. This change in sign is almost 
immediate for QE conditions because p' ~ q, while above the QEP p' > q and for p small 
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now 9^0° and again the response is negative and changes its sign as p increases. An 
interesting feature of this response is that the effect of the FSI is about the same as that 
in the unpolarized case in the QEP and above it, but below the QEP the FSI produce an 
increase of the response that is significant at low q. In fact, for fixed q (i.e., for each row in 
Fig. 3), we see that there is an inversion of the PWIA and DWIA responses at some point 
below the QEP between regions a and b (first and second columns of Fig. 3). At q = 500 
MeV/c and QE conditions (region lib) the effect of the SOP is somewhat bigger than in the 
unpolarized case. 

5.2.2 T responses 

Unpolarized response. The unpolarized response W 00 y ' shown in Fig. 4 presents fea- 
tures that are similar to those seen in the longitudinal one, with the exception that it is 
rather insensitive to the spin-orbit part of the interaction. As we can see in the figure, 
the solid and dashed lines are very similar. An interesting characteristic can be seen if we 

rp/ j \ 

compare the behaviour of the variations AW 00 m both in DWIA (Tables 2-4) and PWIA. In 
DWIA the total variation — relative to the longitudinal — increases with energy for q = 500 
MeV/c from ~ 100% (region Ha) to ~ 104% (region lie), see table 3, while in PWIA we 
have verified that the reverse occurs: there is a decrease from ~ 101% to ~ 91%. Therefore 
one expects an increase of the T/L ratio with u> due to FSI, whereas in pure PWIA the T/L 
ratio decreases with energy. This means that the FSI is stronger in the L-response than in 
the T-response. As we can see comparing Figs. 1 and 4, this is due to the fact that the 
SOP produces an additional reduction of the L-response, while the T-response appears to 
be insensitive to the SOP. 

Vector response The same behaviour of the small SOP effects in the unpolarized trans- 
verse response is also observed in the polarized T responses. For instance, the very large 
effect observed in the ' response does not happen for ' (not shown). We conclude 

that the L responses are sensitive to the SOP, while the T responses are almost independent 
of it. This may be due to the different spin-dependence of the corresponding electromagnetic 
operator. In the L response the electromagnetic operator in leading order involves the spin- 
independent charge operator, which cannot flip the spin. Therefore the transition matrix 
element of the L response is sensitive to the change of the spin part of the wave function 
induced by the SOP. On the other hand, the transverse current is proportional to the Pauli 
spin matrices in leading order, and therefore it is less dependent on the relative change of the 
I + 5 and I — i components of the partial waves. With regard to size, the "* response 

is about 25% of the unpolarized longitudinal for q = 500 MeV/c and is much bigger for 
q = 700 MeV/c (see Tables 2-4). We have observed that the vector transverse response 
approximately scales with the vector longitudinal one when this is computed without the 
SOP and that the scaling factor is about the same as in the unpolarized case (and this is 
the reason for not showing figures for this response). 
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5.2.3 TL responses 

TL(-\-) 

Unpolarized response. The W 00 shown in Fig. 5 is an interference between L and 
T multipoles and thus some similar behaviour to the responses discussed above is expected, 
although the interference can produce a reduction or enhancement of some of the features. 
In fact, in PWIA the TL response is of first-order in S = r\ sin 0, as we can see from Eq. (121), 
and it is due to the interference between charge and convection, and spin-orbit and magne- 
tization terms. In the QEP, the TL strength is around 50% of the longitudinal (Tables 2-4). 
For q = 300 McV/c the TL response is bigger than the T response. This can be understood 
by looking at the expressions for the responses, Eqs. (78 — 83) of reference [6]: in leading 
order the T response is proportional to the single- nucleon response w T ~ 2J^ = 2tG\ i , 
while the TL response is proportional to 

w TL = 2V2( Pc J c + Pso J m )6 = 2V2 (g% + 2G ^ZE E ^G M ^j S. 

For q = 300 MeV/c we have r ~ k 2 , while for the maximum of the momentum distribution 
at p ~ 150 MeV/c, 9 ~ 90°, we have 6 ~ n and therefore 



w TL y/2 

Thus, for low q the TL response is bigger than the T response, although the latter is of 
zcroth-order in 5, while the TL one is of first-order in the expansion. Of course, the currents 
contain terms which are (q, cj)-dependent, whose magnitude has to be taken into account. 
For high q the factor r is bigger and the T response dominates over the TL one. 

With respect to the FSI effect, it is bigger than the one found in the T and L responses. 
The reduction in all cases is more than 50%. This suggest that the convection and spin-orbit 
electromagnetic matrix elements are more sensitive to the dynamics of the FSI. Also this 
response is more sensitive to the SOP than the longitudinal one. In fact, if the SOP is set to 
zero, the FSI effects are almost the same as in the L and T responses. Therefore it appears 
to be the spin-orbit potential that is responsible for the very significant reduction of this 
response when compared with the other two. 

In order to understand this feature better, we have performed a calculation of the re- 
sponses with the SOP set to zero and without spin-orbit and convection pieces in the current. 
As a result, we have checked numerically that the interference between charge and magne- 
tization goes to zero when the SOP docs so. This happens for all of the polarized TL 
responses. In addition, when the SOP is turned on, the zeroth-ordcr contribution is nonzero 
and negative for q > 500 MeV/c at the QEP, providing the reason why the FSI are so large 
in this response. 

TL( — ) 

Vector responses. There is no counterpart of the W w y ' response in the T and L cases. 
We show this response in Fig. 6. As seen in Eqs. (63,64) the vector responses without 
a tilde contribute only to the W TL structure function, which is nonzero only for out-of- 
plane emission (see Eq. (57)). This vector response presents an oscillatory behaviour, its 
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magnitude being ~ 26% of the unpolarized longitudinal in region lib. From the figure it is 
evident that this response is very sensitive to the interaction, the results with and without 
SOP being completely different in all of the kinematical regions here studied. In fact, the 
depencence of this response on the SOP is much more dramatic than for the case of the 
vector longitudinal (compare Figs. 2 and 6), and therefore such big effects, even changing 
the sign in many cases, could in principle be measured if a separation of this response could 
be performed. 

Quadrupole responses. The responses W^J^^ are nonzero in PWIA and, under that 

XVLC+I L( I ) T(-\-) 

approximation, the relation between W 2M and the corresponding W 2 h and W 2 J A is 

TL( -\-) 

the same as in the unpolarized case. The FSI effects in the W 20 response shown in Fig 7 
are stronger than in the case of W^q and W 7 ^, as can be seen comparing Figs. 3 and 7. 
This is again a consequence of the different sensitivities of these responses to the spin-orbit 
interaction (compare the dashed curves in Figs. 3 and 7). The reason for this behaviour 
is again that the SOP generates a contribution from zero-order terms (interference between 
charge and magnetization) in the current matrix elements, which would be zero if there were 
no SOP. 

T Li -\-) 

The FSI effects in the W 20 ^ response also are very dependent on the choice of kine- 
matics. For q — 500 MeV/c at the QEP (region lib), the FSI reduce the response by about 
60% and below the QEP in region Ha, they produce a change of sign. The SOP importance 
is very large at low q, but small above the QEP in region lie. On the other hand, at the 
QEP the importance of FSI and SOP strongly increases with q: at q = 300 MeV/c the total 
effect is about 10%, while the response is drastically reduced for q = 700 MeV/c. 

5.2.4 T' responses 

Polarized electrons are required for the primed responses X", TV . There are a total of 24 
reduced response functions of this kind which we now discuss. Note that there is no T' 
response in absence of nuclear polarization. In fact, for J — even the T" responses are only 
of the tilde type, and therefore all of them are zero for M = 0. 

Vector responses. The W w response shown in Fig. 8 is nonzero in PWIA and has 
no counterpart in the set of unprimed responses, since there is no vector T response of this 
kind. From Eq. (127) we see that in leading order in S this response is expected to be 
of the same order of magnitude as the T responses and in fact in PWIA this response is 
comparable with the T unpolarized one, and the same happens in DWIA, as we can see by 
inspection of Tables 2-4. The FSI produce a reduction that in general is smaller than that 
found for W 00 y ' ; however we find several exceptions to this statement, the most important 
one being that in region Ilia the FSI only slightly affect this response. In this region the 
only noticeable FSI effect is found for parallel emission (8' ~ 0), although for higher missing 
momentum (p' ~ 200 MeV/c) the response is almost independent of the interaction. In 
addition for q > 500 MeV the SOP effect is negligible. Therefore this could be a good place 
to look for other model dependences such as on the specifics of the single-nucleon current or 
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nuclear structure. A second exception is in region la, where the FSI enhance and produce 
a shift of the response, although we are far away from the range of applicability of the 
potential. 

With respect to the W n response (Fig. 9), one would expect that it should in general 
be bigger than W n y ' which is zero in PWIA, but we find that in many cases the FSI produce 
a significant reduction of the PWIA results for this response, the effect being more than 60% 
in the QEP region, where the variation AVF X1 of this response is always smaller than the 
corresponding AW n (see Tables 2-4). On the other hand, in regions IIa,c, this variation 
is <~ 50% of the unpolarized longitudinal case, and it is above 100% in regions IIIa,c. The 
reason why this response is so small at the QEP can be seen by looking at its behaviour 
for fixed q, as a function of uo (Fig. 9). There is a tendency for this response to change 
from positive at lower ui to negative at some point slightly above the QEP. Comparing the 
DWIA and PWIA results we also see a general tendency for the former to be reduced in 
magnitude with respect to the latter, where this reduction is seen to be significant below 
the QEP and is larger for low q. On the other hand, above the QEP (regions Ic, lie) the 
mean (negative) peak of this response is almost not affected by FSI, so these could be good 
places to investigate other aspects of the reaction. This trend does not apply for q = 700 
MeV/c, uj = 300 MeV (region IIIc), but here again we are far away from the validity of the 
potential parametrization. As in the case of the T response, this response is insensitive to 
the SOP. 

T' ( — ) T'( — \ 

Octupole responses. The W 30 response shown in Fig. 10 is comparable to the W w 
one in both PWIA and DWIA, and over the entire kinematical range considered. In some 
places there are drastic changes due to the FSI with respect to the PWIA. For example, in 

T' ( — ) 

region Ha below the QEP this response is almost opposite to the W 1Q in PWIA, that is 

reduced by FSI by less than 20%, while the W 3CI is almost canceled by FSI. On the other 
hand, in region lib, the FSI effect provides a reduction of around 20%. If we go further in 
energy (to region lie), where the PWIA presents an oscillatory behaviour, the FSI yield a 
reduction by 50% at the maximum, and almost do not affect the first minimum. We observe 
a similar behaviour in region Ic where the reduction due to FSI is even bigger, around 75%. 
Again, the SOP effect is negligible in this response. 

The W 31 response is shown in Fig. 11 and is also important. It reaches its greatest 
values above the QEP (see also Tables 2-4). For instance, its variation is about 25%, 50% 
and 120% of the unpolarized longitudinal in regions Ic, lie, and IIIc, respectively. At the 
QEP, its magnitude is around 70% of W n , and the FSI effects are almost identical in both 
of them, in general very strong for low q (producing even a change of sign) and moderate 
at q = 700 MeV/c (with a reduction of ~ 65%). Below the QEP in region Ha this response 

is around 80% of the W lx and about one half of W 10 , the FSI effects being smaller 

than in the W u response. FSI effects are also seen to be very different at the QEP 
upon comparing the M. = and A4 = 1 odd X" responses. It is important to note that the 
identification of places where FSI effects are very different, as in the case of the W n y ' and 
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W 31 in region Ha, could play an important role in this type of study where the goal is 
to isolate such effects. The remaining octupole responses W 32 and W 33 1 ' (not shown 
in the figures) tend to be smaller than the one discussed above, although the former could 
be noticeable in some cases, while the latter is negligible. 

In all cases we find that the influence of the SOP in the T" responses is small, as in the 
T case, thus making it possible to use these transverse responses to study the central part 
of the interaction and to use the responses which are sensitive to the SOP, such as the L or 
TL cases, to study the spin-dependent part of the interaction. 



5.2.5 TL' responses 

T ' L' ( I ) 

Unpolarized response. The "fifth response function" W 00 (shown in Fig. 12) is zero 
in PWIA. It can be separated by performing an asymmetry measurement by flipping the 
electron helicity in experiments with unpolarized nuclei, therefore providing an alternative 
way of studying the FSI. First we note that its size is large and comparable with the un- 
polarized TL response, with which it is related by Eqs. (46-49). In fact, the TL and TL' 
responses differ in that they select different combinations of the real or imaginary parts 
of the interference multipoles in Eqs. (28,29). Specifically, the W 00 1 ' and W m 1 ' con- 
tributions contain the real and imaginary parts, respectively, of the same combination of 
multipoles (see Eqs. (37,38)). Therefore their magnitudes are roughly determined by the av- 
erages (cos AS) and (sin AS) of cos[Re((5;/j' —Sij)} and sin[Re(<Jj'j' —Sij)], respectively, that is, 
of the differences between the real phase-shifts for different partial waves. The magnitudes 
of these averages are determined by factors involving the weighting amplitudes IC^iE^l and 
\C a /M a \ in Eqs. (46-49). The intercomparison between TL and TL' responses can therefore 
provide information about the mean real phase-shift deviation. Of course, within the aver- 
age the relation (cos AS) 2 + (sin AS) 2 = 1 is in general not true, although for QE conditions 
it is approximately satisfied. In fact, the values of the pair (Wo L{+ \Wo L ' {+) ) at the peak 
in regions lb, lib and Illb are - (70,40), (70,45) and (30,70) [GcV^ 1 x 10~ 3 ] and therefore 

r tl(4-)~\ r TL'(+)~\ 2 

the combination W Q0 + W QQ ~ constant. Outside the QEP this relation does 

not hold, but in any case the two responses are always comparable in magnitude. Moreover, 
the effect of the SOP is of similar importance, ~ 20% in both responses. 

With regard to their behaviour as functions of the kinematic variables, the TL and TL' 
responses have different structures, the most significant variations being found below the 
QEP and at low energy. It is interesting to note that the behaviour of the W 00 1 ' response 

— L(— ) — 'T(-) 

is more similar to that of the and W n responses. Our results show that the effect 

of the SOP in the W 00 and W n responses is similar, and it presents a tendency to 
modify the response in the same way. It is also noteworthy that when the SOP is set to zero, 
the three responses Wq Q L W^" 1 and acquire the same shape and that W^ L ^ 

also shows the same tendency (compare Figs. 2, 6 and 12). This is related to the fact that 
all of these responses are given by imaginary parts of the single transition responses for each 
partial wave. 
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Vector responses The TV vector responses are large in magnitude and in all cases 
considered they are comparable to the Wqq + ^ response. From Tables 2-4 we see that their 
importance increases with the momentum transfer, and for high q the W^[ L response is the 
dominant one, being around 160% and 220% of the unpolarized longitudinal case for q = 500 
and 700 MeV/c, respectively. The W w response is shown in Fig. 13, and it presents a 
behaviour which is opposite in sign to W n y , both with and without FSI, but the effect 
of the SOP is different at the QEP (regions lb, lib, Illb). For instance, in region lib the 
introduction of the SOP changes the sign and peak position, and in region Illb the response 
is negligible without SOP, while its inclusion produces an increase of the absolute strength 

of the PWIA results by more than 100%. On the contrary, in region Ic the FSI effects are 

TjJ i ^ — ~~tl'( ) 

negligible. The W n and W w responses are very large in magnitude and almost 

— tl' (—) 

opposite in sign. It is interesting that the W u response, not shown in the figures, has 
a behaviour that is very similar to W^ + \ although opposite in sign, and that the effect of 
the FSI and SOP is also similar in these two responses. This is related to the fact that all of 
these responses are composed of real parts of responses for single multipolc transitions R^, a , 
R T J a , R T J- a , R^ 1 , and R T J£ in Eqs. (40,42,46,48), and that all of them are dominated by 
the zeroth-order static component of the electromagnetic current. 

Quadrupole responses. The W 20 , shown in Fig. 14, has no T' counterpart. This 
response is bigger than the fifth response function below the QEP, but is smaller at and near 
the QEP. This response shows a large sensitivity to the SOP, which completely changes its 
behaviour in many cases, especially at the QEP. Specifically, the SOP produces the tendency 
to reduce this response to negative values for quasielastic conditions. In general all of the 

TL' f-f-l ~~^ r TL r (-\-) 

W 2M y ' and W 2M y ' responses have almost the same magnitude in region lib, and they 

are about twice that of the transverse W 2 A4 ( see Tables 2-4), which could help in making 
an experimental separation of these quadrupole responses feasible. In general all of the 
quadrupole TV responses are very sensitive to the SOP, in contrast to what happens for 
the quadrupole T" responses. 

Octupole responses. There arc four W^J ^ and three W 3 j^ ^ octupole responses and 
in general all of them are big compared with the fifth response function, the biggest ones 

TL 1 ( — ) 

being for M. = 0, 1, as seen by inspection of Tables 2-4. In fact, W 30 , shown in Fig. 15, is 

T'( — ) 

the response with the biggest variation in regions IIa,c and IIIa,c, and W 30 is the second 

response of importance in this ranking. The W 30 L W 31 L ^ and W 32 L ^ responses 
are strongly affected by FSI, both in size and sign at QE conditions, and the last two are 
significantly suppressed by the FSI. The SOP does not play a particularly important role in 
these responses, with some exceptions. In regions Ha, lie and Ilia, the W 30 y ' response is 
larger in magnitude than the unpolarized longitudinal response even with the inclusion of 
FSI and it is rather insensitive to the SOP. Also W 31 is large in these regions, where 

T'(-) 

it is similar in shape and magnitude to W 30 . In summary, the large magnitude of these 
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octupole responses indicates that they could possibly be separated through a sequence of 
asymmetry measurements. 

5.3 Relativistic corrections 

In this section we briefly discuss the importance of the different terms in the first-order 
expansion of the relativistic current in powers of S = r\ sin 9 [6]. In particular, the goal is to 
study the importance of the first-order terms which we call spin-orbit (in the longitudinal) 
and convection (in the transverse), but which differ from the terms usually used in non- 
relativistic calculations in that they maintain the full dependence on (q, oj) in the relativistic 
current (sec the definitions in Eqs. (123-126)), and therefore they are more suited for use at 
high momentum transfer near the QEP. Also the spin-orbit charge (SOC) term is usually 
not included in non-relativistic calculations, but as seen in Eq. (121), it contributes to the 
single-nucleon TL response at the same level as the convection current and therefore its 
inclusion in the longitudinal component can be important for properly describing the TL 
responses at high momentum transfer, as was found in Ref. [6] in the inclusive polarized 
case. Moreover, the FSI can affect each one of these contributions in a different way, as may 
be reflected especially in the responses which are zero in PWIA. 

In the following we show only two figures to illustrate the most interesting cases. As in 
the last section, we show a separate response function in the nine panels corresponding to 
the same choices of kinematics made above. Now the solid line in Figs. 16-17 represents the 
total response computed in DWIA with the complete optical potential. With dashed lines 
we show the responses computed in leading order in the expansion, and therefore the dashed 
lines only include charge plus magnetization (zeroth-order) terms. Finally, the dot-dashed 
lines include, in addition, the convection current, and thus they are representative of what 
can be found in a standard non-relativistic calculation (with the important difference that 
we include the relativistic (q, u>) behaviour which is significant at high q). In other words, 
the effect of the spin-orbit term in the longitudinal current is to transform the dot-dashed 
curves to the solid ones. 

We now discuss the importance of the different parts of the electromagnetic current in the 
case of the TL responses. In PWIA these responses are entirely due to interferences between 
charge and convection, and spin-orbit and magnetization terms, and as a consequence the 
first-order terms of the current are crucial when describing these responses. Of course, when 
FSI are turned on, the cancellation of the zeroth-order terms is not exact and the interference 
between charge and magnetization is nonzero with FSI (dashed lines in the figures), although 
its magnitude is comparable with the total response or smaller in some cases. As stated 
above, we have checked that the zero-order contribution to these responses goes away when 
the SOP is set to zero, and this is the reason why the TL responses are so sensitive to the 
spin-orbit part of the interaction. 

The W 00 y , shown in Fig. 16, is especially sensitive to the first-order terms. The 
zeroth-order terms (interference between charge and magnetization pieces) are shown with 
dashed lines and they are zero in the absence of a spin-orbit interaction in the optical 
potential. The importance of the zeroth-order terms increases with the momentum transfer 
in all kinematic regions. For q = 300 MeV/c this contribution is almost negligible, while 
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for q = 700 MeV/c it is of the same order of magnitude as the total responses and opposite 
in sign. The increasing importance of the zeroth-order term is closely related with what 
we found in last section — the effect of the SOP increases with q in this response and the 
inclusion of the SOP potential in the FSI prevents the spin cancellations from occurring in 
the zeroth-order term. Shown with dot-dashed lines is the response computed including only 
charge, magnetization and convection terms in the electromagnetic current. These are the 
expected results if a standard non-relativistic treatment is used to compute the responses 
(no spin-orbit correction to the longitudinal piece) . The effect of the convection term is to 
completely change the behaviour of this response with respect with the zeroth-order, even 
changing the sign for q > 500 MeV/c. When we include also the SOC piece, we obtain the 
solid lines in the figure, thereby showing the importance that this correction has in the TL 
response. The SOC contribution also increases with q. The effect of this term is to increase 
the total (dot-dashed) response by ~ 15%, <~ 75% and <~ 200% at the peak for q = 300, 
500 and 700 MeV/c, respectively, and QE conditions (regions lb, lib, Illb). The effect is 
even bigger below the QEP, in regions Ia,b,c. As seen in Eqs. (124,126), the spin-orbit 
charge is proportional to nj \J\ + r, which increases with q, while the convection piece is 
proportional to y/r/ k, which slightly decreases with q. Therefore the relative importance of 
the SOC versus the convection current increases with q, as seen in the plots. Accordingly 
we conclude that any non-relativistic model which does not include the spin-orbit term 
significantly underestimates this response at a typical value of q = 500 MeV/c and QE 
conditions. 

The zeroth-order term is important in the W 1M responses. For QE conditions, the 
convection effect is crucial at q = 300 MeV/c, and decreases with q. In region Illb spin-orbit 

and convection effects are only around 20% of the total response in W 10 and smaller in 

TL ( — ) 

W n y ' . These effects are much more important below the QEP, although in this region 
these responses are small compared with the Wqq + \ On the other hand, the first-order 

— 'TL(-) 

terms are more important in the W n case shown in Fig. 17, which is exactly zero in 

L( — ) 

PWIA. The zeroth-order terms are small in W n which presents the same trends as the 
unpolarized response, Wq L ^ + \ and both responses are of the same order for q = 700 MeV/c 
(see Table 4). For instance, in region Illb at the QEP, these two responses are the most 
important ones of the TL-type. The spin-orbit effects are 20%, 60% and 75% of the response 
computed with only zeroth-order plus convection terms in QE conditions at the peak. 

Similar trends as those found in the unpolarized TL case are observed for the rest of the 
TL-responses although each one of them presents particular features for different kinematics 
and different values of J M. . For obvious reasons it is not possible to show all of the figures 
here; the complete set is presently available from the authors at [21]. 

We now discuss the TT responses very briefly because they are also in general very small 
and hard to measure. The unpolarized response is mainly due to the zeroth-order magne- 
tization term, the convection current contributing only a small fraction to this response. 
The PWIA response is completely due to the convection current and the response com- 
puted in PWIA is much larger than the effect of the convection upon the DWIA response, 
which implies that the pure convection current is significantly suppressed by the FSI. Also 
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the response without the SOP is more similar to a pure reduction of the convection PWIA 
results than with it. This suggest that the SOP significantly modifies the cancellation of 
the pure magnetization response. And in fact, when a calculation was made of the different 
contributions to the TT responses with the SOP set to zero we found that the magnetization 
contribution goes away for all of the TT responses, and not only for the unpolarized one. 
In the other responses of the TT-type, we find that the magnetization contribution is in 
general dominant; yet in many cases, the convection contribution is important and several 
times crucial to describe these responses properly. The effect of the convection current in 
general decreases with q. 

More interesting is the effect of relativistic corrections upon the fifth response function. 
In general the first-order terms (spin-orbit and convection contributions) are found to be 
negligible in the QEP region and above it, although there is a noticeable effect below the 
QEP in regions la and Ha where the fifth response function has two maxima at p ~ 80 and 
220 MeV/c, and one (positive) minimum in between at p <~ 130 MeV/c (see Fig. 12). The 
effect of the convection current is to decrease the response for low p and to increase it for 
high p. The spin-orbit charge increases this response for low p, while it slightly decreases 
it for high p. The effects are noticeable: for instance, in region Ha the two maxima would 
be approximately equal if only magnetization is present. If the convection also contributes, 
then the second peak is bigger than the first one by a factor 1.5. If the SOC is also present, 
then the difference between the first and second maxima is about 20%; that is the situation 
shown in Fig. 12, where the complete current is used. Of course, the FSI effects play an 
important role here and the situation is completely changed if the SOP is zero. 

Finally, we have found non-negligible first-order contributions to some of the remaining 
TL' responses, the most important case being the W w 1 ' and W 20 y ' responses for a 
variety of different kinematics. 

6 Conclusions 

In this work we have investigated polarization observables in coincidence electron scattering 
from nuclei when both beam and target are polarized. Our primary objective has been to 
develop the formalism for the A(e, e'N)B reaction, introducing a complete set of reduced 
response functions which are independent of the polarization angles and azimuthal nucleon 
emission angle, and which contain the maximum information about nuclear structure and 
reaction mechanisms of interest. We have applied the above formalism to the particular 
case of a one-hole nucleus in the shell model using both PWIA and DWIA, and we have 
examined theoretically the differences which arise when the nucleon is ejected from different 
complete or incomplete shells, and when the daughter nucleus is left in states with different 
angular momenta. With respect to the theoretical results for the extreme shell model it is 
worth pointing out the following: 

1. The reduced response functions for nucleon knockout from the same shell but different 
spin for the daughter nucleus are proportional to the responses of a single polarized 
particle in a shell. The proportionality constant contains a recoupling coefficient of 
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the angular momenta involved. When a more sophisticated nuclear model with con- 
figuration mixing is used, interferences between different shells are expected [7] and 
deviations from this proportionality could provide information about the nuclear struc- 
ture involved. 

2. Nucleon knockout from a complete inner shell in a polarized nucleus can show po- 
larization effects although the involved shell has spin zero. The reason is that the 
final nuclear state is left in a state with definite angular momentum, and therefore the 
remaining incomplete shell is partially polarized. When the sum over final daughter 
spin is performed, all of the polarization responses add to zero. Again, deviations from 
this result arise in a more realistic model and, for instance, for a deformed nucleus 
such as 21 Ne, it has been found in Ref. [7] using PWIA for the ejected proton that as 
more states of the residual nucleus are involved, the effect of the target polarization is 
weakened. 

3. The two above results are also valid for PWIA with the addition that similar con- 
clusions are obtained for the polarized momentum distribution of the nucleus: it is 
proportional to the polarized momentum distribution of a single hole in the shell from 
which the nucleon is ejected, with the same proportionality coefficient as found be- 
tween the reduced response functions. Within the context of this approximation, the 
study of this reaction opens the possibility of measuring the complete tri-dimensional 
momentum distribution of the nucleus. 

In this work we have focused on the study of the angular reduced reponse functions 
because they are the basic quantities that describe the reaction. Nevertheless, the large 
number of these functions (72 for the case of 39 K studied here) makes it very difficult to 
attempt a separation of all of them experimentally, although, of course, specific combina- 
tions of the reduced responses could be isolated by performing a combination of asymmetry 
measurements. To understand the relative importance of these combinations in (the usual) 
circumstances where they are not entirely isolated it is important to perform an analysis 
of all of the response functions and to study their sensitivities to details of the reaction 
mechanism prior to any measurement. 

That is what we have done in this work for the particular case of the 39 K nucleus described 
as a hole in the d 3 / 2 shell of 40 Ca. For the ejected proton we have used the DWIA with 
an optical potential and compared the results obtained with the PWIA. We have included 
relativity in the kinematics and in the electromagnetic current. 

We have studied the model dependence of the reduced response functions by intercom- 
parison between models, specifically addressing the question of what is the role of the central 
and spin-orbit pieces of the optical potential and exploring the effects of the first-order pieces 
of the electromagnetic current (convection and spin-orbit charge) in the responses. Let us 
summarize the main conclusions that may be drawn from our work. 

1. In general the FSI produce a reduction of the absolute strength of the unpolarized 
responses with respect to the PWIA. The percentage of this reduction depends on 
the kinematic conditions, being typically of <~ 35% for the L and T responses and 
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for q = 500 MeV/c at QE conditions, however being more than a 50% effect for the 
TL response. The reason of this big reduction is because the FSI generate a nonzero 
contribution from the charge-magnetization, zeroth-order interference term, which is 
exactly zero in PWIA, but is negative with FSI for q > 500 MeV/c. We have checked 
that this effect is produced mainly by the spin-orbit potential, since the zeroth-order 
terms go away when we set the SOP equal to zero. The TT response is a few percent 
of the others and negative in PWIA, but again the FSI produce a drastic change, even 
inverting the sign, mainly due to the zeroth-order magnetization term in the current, 
which is zero in PWIA, although not in DWIA when the SOP is turned on. 

2. The fifth response function Wq L ^ which arises for unpolarized nuclei and polarized 

electrons is comparable in size with the Wq L ^ response, being around 20-90% of the 
longitudinal unpolarized response under different kinematics, and it is quite sensitive 
to the SOP, the most important effects being seen below the QEP. 

3. In general the FSI effects are large in those responses that are nonzero in PWIA, 
especially the quadrupole L, T and TL, and vector and octupole T 1 and TL' cases; 
there is dependence on the particular rank J M. of each response and in some cases 
the FSI even produce a change of sign. The variety of effects found is too large to be 
summarized in a few lines, but we can cite as most outstanding cases with dramatic 
dependence on FSI the responses W^ +) , W^ +) , W^ L{+) , w£ L(+) , W^ ( '\ W^~ ] , 

TL' ( ) T L' ( ) 

Wio y ~' and y ~>. Also some cases have been found where the FSI are negligible, 
therefore providing an ideal place to study other reaction issues, such as those involving 

T'( — ) 

the nucleon current or details of nuclear structure. The octupole responses Wg .]j 

and W 3M are found to be large for M. = 0, 1, especially for high q, and these 
responses are dominant over the others below and above the QEP, probably making 
their extraction experimentally feasible. 

4. The polarized responses which are zero in PWIA are a fortiori sensitive to the FSI. 

— L(-) — T( — ) 

Additionally, the W n response is rather dependent on the SOP, although the W n 
response is not. The effects of the SOP arc dramatic in W 1M y ' due to the fact that 
zeroth-order terms in the current only contribute to this response if the SOP is nonzero. 

TL' ( I ) 

Also the W 20 y ' is very sensitive to the SOP. 

5. All of the T and T' responses are largely insensitive to the spin-orbit interaction; 
therefore they could be used to study the central part of the FSI. On the other hand, 
the L, TL and TL' responses are much more sensitive to the SOP, especially in the 
case of the TL responses. 

6. The spin-orbit first-order term in the charge operator is crucial for properly describing 
the TL responses and some of the TL' ones, especially at high q where relativistic 
effects are expected to be more important. At q = 700 MeV/c the spin-orbit and 
convection contribution are of the same order for both polarized and unpolarized TL 
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responses. In other cases the first-order effects are smaller, but noticeable in some 
situations, for instance, in the fifth response function. 



7. The non-locality correction with a Percy-like factor is in general smaller than the 
effects mentioned above and it can be safely ignored in a first exploratory study of this 
kind. 

In conclusion, we have found a wide variety of effects of the FSI and relativistic corrections 
in the different reduced response functions of polarized nuclei. As the effects are different for 
different responses, systematic intercomparisons between selected separated responses could 
provide important information about different aspects of the reaction mechanism. The 
quantity of information which may be obtained with exclusive polarized electron scattering 
from polarized nuclei is more than one order of magnitude richer than in the unpolarized 
case and this should allow much more stringent constraints to be placed on any particular 
model of the dynamics. 

Experimentally, however, the amount of work involved in measuring and extracting all 
of the angular reduced response functions is daunting, with little hope for a complete ex- 
traction of all of them in the near future, and therefore it is desirable to extend this analysis 
to the study of spin observables which can be directly measured, such as cross sections and 
asymmetries. Our aim here has been to provide a first detailed insight into the basic quan- 
tities describing the reaction, namely the reduced response functions, as a first step towards 
understanding the role of the different reaction ingredients in the problem. It is our intent 
to apply the present formalism to study various classes of asymmetries and investigate their 
utility for extracting fundamental properties of polarized nuclei such as polarized momentum 
distributions and spin distributions. Work along these lines is in progress. 
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Appendix A. Calculation of the longitudinal response 

Here we present a proof of Eqs. (34,40,41) for the longitudinal response. We begin from the 
definition, Eq. (9), and inserting the multipole expansion in Eq. (15), we have 




(132) 



J J' 
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where Bj)j is given by Eq. (17), with T'j, = Mj> and Tj — Mj. Using the general result in 
Eq. (21) we obtain 

K L =**T, E i J - J 'i l '- l ^v + j-fj\rjYj'U*°>A o "J' ) (133) 

<t<t' JJ'I ^ ' 

where the Coulomb amplitudes in Eq. (22) have been introduced. Now we use the fact 
that the Coulomb operator Mj has natural parity (— 1) J . If we denote by IT and Hb 
the parities of the target and daughter nuclei, respectively, we conclude that the matrix 
element C a is nonzero only for transitions satisfying IIjIlB(— l) l+J = 1, while for C„i we 
have n 4 n B (-l)''+ J ' = 1. Therefore, for the product C*,C a we have the selection rule 
I + J + I' + J' = even. Taking this into account and using the parity function in Eq. (33), 
we can write 

iJ-i-V-nCZC = ti^j^CZC (134) 
and we have for the longitudinal response 

* L = 4*E E Pi + + v + j'fj[yjYj'\LO^J o o L A$.-v,j-iC**.C a . (135) 

era' JJ'L ^ ' 

Now, using £+„ = £+ m , we can write the sum over a' a in a symmetrized way 

E' P W'+J' $<J ' ,T ( ) €j'-l',J-l C <r' C(7 = 
<j<j' ^ ' 

\Y. P ^+j\i o )^-v,j-i[^C:'Ca + ^C*M, (136) 

because, due to the parity functions, I + I' + J 1 = even, I + I' + J + J' = even, and due to 
the parity properties of the 3-j, we have that J + J' + L = even. Using now the symmetry 
property in Eq. (20), we have $ CT(T ' = {—l)^$ a i a because, due to the parity functions, we 
have that J' + L = even. Therefore we can write the symmetrized product as 

\ [$ ffV c;,c a + ^c:c a ,\ = U> a , a [c;,c„ + (-i) J c;c^] 

= ^[P+R^ + iPjlX], (137) 

where we have used the definition in Eq. (25), and noticed that when J = even, we are 
taking twice the real part of C*,C a , and when J = odd, we have twice the imaginary part 
of C*,C a times the imaginary unit. 

On the other hand, we have the conjugation property 

[YjYj,}* L0 = {-lf +J ' +L [YjYj,U = (-1) J [YjYj,} L0 (138) 

because again J' + L = even, and therefore the coupling [YjYj>]l is real for J = even and 
imaginary for J = odd; and so we can write 

[YjYj,] L0 = P+Ajj, lq + iP-Bjj, L0 . (139) 



41 



Finally, inserting Eqs. (136,137,139) in Eq. (135) we obtain: 

k l = ^EEW^f! o o)ow* 

a 1 a JJ'L V / 

x [P+Ajj, L0 RL, a - PjBjj, LO lL, a ] , (140) 

from which it is straightforward to read off Eqs. (34,40,41). Although the derivations of the 
remaining responses are somewhat more involved, the basic steps for obtaining them are 
natural extensions of the ones given here. 



Appendix B. Reduced responses for a one- hole nucleus 

In this appendix we show how to perform the sums over the angular momenta of the final 
states in the reduced responses for a one-hole nucleus whose ground state can be described 
as an extreme one-hole configuration Since the following steps are analogous for each 

one of the reduced responses Wjj^ l , we only consider the case of the Wj [ p L response, 
given in general by Eq. (40) 

W L j { J>l = Y. P ^+J>^'°{1 oW-W^ (141) 

era' ^ ' 

with R^, a = ReC*/C<j. First consider the case h ^ i, when the struck nucleon is in a shell 
other than the polarized one. In this case the final states are 

. \C), (142) 

where \j) = \E'lj) is a single particle in the continuum with energy E' . Here it is useful to 
define, in addition to the index a = {l,j,jf,J}, indices 7 = {l,j,J}, 7' = {l',j',J'}, that 
do not contain the final angular momenta jf, j'j which are going to be summed. With these 
indices, the many-body Coulomb matrix element can be written as 

= (fWMAr 1 ) = [JB]\jf](-l) J - j - jh { j j j / h j B } c 7 , (143) 

where we denote with a lower-case letter the Coulomb matrix element for the single-nucleon 
transition from the shell h to the continuum: 



l/)= 4[b{b\}j B 



c 7 = {j\\Mj\\h). (144) 
We can then write the real part of the product C*,C a in the following way 

(145) 
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where now Ty 7 is the real part of the interference between single-particle Coulomb multipolcs 

r£ 7 = Rec*,c 7 . (146) 
Therefore we can write the response as 

jfj' f i'i 

x(-l)^'W+i( 3i f f J' || ji if J \ r L (14?) 
I 3 Jh Jb J i 3 Jh Jb J 77 

Now we are interested in performing the sums over jf,j'f- Specifically we are going to 
compute the following coefficient: 



Inserting the definition in Eq. (19) of <fr a i ai 

jfff 

f 

_L 

2 

x "«'i? 2 i Hj a i| (149) 

The sum over jf 7 j'j can be written as the product of a 6-j and a 9-j symbol using the 
relation 




a' V C c'\\ d h 6 f \ 
d e f J I h 1 3 J 



— i ■ A I j /i A a I 1 /i 7 a . . 



(150) 



with 5 = & + c — a' — d — d — e + 2h + j. Note that Rotenberg's version [13] of this 
equation (Eq. (3.23)) has an error in the third 6-j of the righthand side: instead of h, in 
that reference appear an i. In order to obtain the correct expression, in this work we have 
proven the above relation using the graphical method of Ref. [14]; for brevity we do not 
reproduce the procedures used here. 
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Using the above relation and after some rearrangement and simplication of phases we 
obtain for the sum A the result 

a = [j}[j'm\j'}[j'}[L}(-^ J ' +jB+J+J ' ( ( \ J q ) 

x(-i) J + J '+^'(-i)^+«-vJ I J jh J : f j.(i5i) 

! 3 f J'\ 

The 9-j coefficient has the same structure as a polarized particle with spin jh would have. 
In order to relate the result with the response of a polarized particle in the shell h and 
spin-zero daughter nucleus, Jb = 0, we first notice that the coefficient $ 7 / 7 corresponding 
to Eq. (19) for that process can be written 

^(jj'L,h^o) = iJ}[JW}[j'mm'}(-i) J ' +j+1/2+J+J ' 

J J' L 

3h 3h J }■ (152) 



f 3 J' \ 3' 3 J' 

1 _ l 

2 2 



1 o J\ j j' 



Inserting the simple result (— l)- 7-1- -? +J for the 6-j symbol with one argument equal 

to zero, the coefficient 4> 7 / 7 becomes 



Sy 7 (.7A^0) = [J][J'][m[J'}[L](-^ J ' +i+1/2+J+J ' 



r 3 j' 

1 - 1 



2 2 



Then we obtain for the sum A the value 



J 


J 1 


L 


3h 


3h 


J 


3 


3' 


J' 


k 


J 




3h 


Jb 


}■ 



(153) 



(154) 



a = {-i) J + J ' +J +i + i'(-iy h+ i i+jB+1 %, 7 (h -» o) ( j : 

[ 3h 

Finally, inserting this value into Eq. (147) we have for the longitudinal response 

W L(+) = { _ l)jh+ji +J B +J [J s] 2 j ji 3i J 

JJ L I 3h 3h Jb 

x E Pl+v+j'**^ ~+ °) ( o o o ) &-i',J-rii ( 155 ) 

7'7 

from which we directly obtain Eq. (76). 

For the other case, h = i, the final hadronic states are of the type 

|/) = -L \a] [blb\} jB ] \C), (156) 
V2 L j j/ 
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and the transition matrix element is now 



C <T = (/||M J ||r 1 )=^[J s ]b'/](-l) J -^{ * | j^}c 7 . (157) 

The rest of the derivation is completely analogous, with the exception of the factor y/2, that 
becomes a factor of two in the response. Therefore, making h = i in the above equations 
and multiplying by two, it is straightforward to obtain Eq. (81). 
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w[MeV] 


J 


M 


W L 


W T 


W TL 


W TL 


W TT 


W TT 


W T ' 


W TL ' 


W TL ' 


40 








100.0 


38.9 


53.3 




10.8 






23.8 






1 









30.5 




7.4 




42.0 


16.6 






1 


1 


29.6 


15.3 


17.9 


34.1 


8.3 


10.5 


11.3 


116.2 


78.7 




2 





129.1 


36.0 


29.7 




9.4 






61.6 






2 


1 


24.0 


4.0 


24.1 


7.7 


5.8 


9.6 


3.4 


48.7 


74.3 




2 


2 


13.7 


7.6 


8.6 


10.1 


3.0 


4.2 


8.1 


4.9 


11.7 




3 









15.4 




8.8 




28.1 


52.7 






3 


1 


15.4 


2.1 


17.1 


13.4 


2.7 


1.2 


10.0 


38.1 


39.2 




3 


2 


9.8 


3.1 


4.7 


5.7 


3.6 


2.2 


2.1 


7.0 


5.7 




3 


3 


2.4 


0.3 


2.3 


1.7 


1.0 


1.1 


1.2 


5.8 


5.4 


56 








100.0 


39.8 


53.6 




5.2 






58.2 






1 









24.3 




3.6 




41.5 


36.6 






1 


1 


55.6 


22.2 


12.2 


31.7 


9.1 


8.1 


10.4 


79.3 


70.7 




2 





89.1 


23.5 


36.5 




9.2 






14.1 






2 


1 


19.7 


6.9 


6.8 


3.3 


4.7 


5.7 


7.2 


24.7 


27.8 




2 


2 


18.9 


8.1 


10.8 


3.9 


1.3 


1.3 


3.9 


14.7 


21.4 




3 









17.0 




7.5 




20.5 


69.5 






3 


1 


7.3 


2.9 


12.4 


9.9 


3.6 


5.0 


10.0 


34.1 


21.5 




3 


2 


7.1 


1.8 


5.0 


3.3 


1.4 


1.8 


4.0 


8.9 


5.0 




3 


3 


0.9 


0.3 


2.2 


1.9 


1.0 


1.1 


0.9 


6.4 


6.3 


80 








100.0 


36.8 


45.1 




3.2 






38.4 






1 









12.4 




3.2 




38.6 


86.0 






1 


1 


38.9 


14.6 


13.3 


21.0 


5.2 


3.5 


26.2 


97.0 


93.2 




2 





99.6 


31.9 


28.3 




2.3 






20.2 






2 


1 


49.5 


17.0 


19.7 


0.9 


1.0 


1.0 


5.3 


22.6 


18.5 




2 


2 


17.1 


6.1 


9.6 


0.9 


1.1 


0.5 


2.0 


11.6 


11.6 




3 









40.3 




15.9 




55.5 


147.9 






3 


1 


4.4 


0.8 


23.5 


14.6 


5.3 


7.0 


23.6 


87.1 


41.5 




3 


2 


0.3 


0.3 


5.0 


3.0 


1.5 


3.0 


5.2 


16.6 


12.6 




3 


3 


0.9 


0.3 


1.5 


1.5 


0.5 


0.8 


1.2 


3.5 


4.6 



Table 2: Total variation (maximum minus minimum) of the reduced response functions as 
a percentage of the variation of the unpolarized longitudinal. The responses are computed in 
D WIA for q = 300 MeV/c. 



47 



w[MeV] 


J 


M 


W L 


W T 


W TL 


W TL 


W TT 


W TT 


W T ' 


W TL ' 


w TL ' 


110.0 








100.0 


100.1 


43.9 




9.8 






20.8 






1 









16.5 




16.1 




74.1 


89.3 






1 


1 


16.0 


25.7 


14.1 


14.1 


16.2 


5.5 


46.6 


70.0 


165.6 




2 





36.9 


44.1 


19.0 




4.8 






51.8 






2 


1 


25.0 


30.8 


21.0 


5.6 


3.0 


4.2 


6.4 


22.3 


28.3 




2 


2 


22.8 


21.2 


9.7 


2.0 


2.9 


5.2 


7.0 


5.5 


7.3 




3 









16.0 




12.3 




25.6 


135.4 






3 


1 


5.5 


5.5 


13.2 


7.6 


3.1 


4.4 


37.7 


35.8 


17.2 




3 


2 


0.7 


2.7 


1.6 


3.7 


2.5 


2.4 


3.9 


23.4 


12.7 




3 


3 


1.5 


0.3 


1.3 


1.5 


0.8 


1.1 


3.7 


7.9 


9.6 


133.5 








100.0 


103.1 


48.7 




9.1 






37.5 






1 









26.1 




8.5 




96.4 


11.6 






1 


1 


23.9 


32.0 


4.9 


17.9 


9.3 


8.0 


18.2 


97.2 


158.2 




2 





50.4 


57.5 


18.8 




9.1 






19.9 






2 


1 


8.6 


8.6 


13.3 


3.9 


3.4 


5.5 


7.4 


19.6 


17.7 




2 


2 


23.7 


24.4 


12.2 


1.4 


2.6 


3.5 


5.3 


13.4 


14.3 




3 









14.1 




5.5 




61.8 


42.9 






3 


1 


1.8 


1.6 


4.7 


2.6 


1.9 


3.5 


17.6 


67.6 


22.9 




3 


2 


2.2 


2.1 


2.6 


3.2 


1.3 


0.9 


10.4 


6.1 


3.7 




3 


3 


1.0 


0.7 


0.8 


0.6 


1.1 


1.0 


1.8 


11.3 


11.1 


160.0 








100.0 


104.0 


50.5 




4.8 






42.9 






1 









22.5 




3.2 




100.3 


95.4 






1 


1 


29.2 


33.8 


12.4 


19.5 


8.0 


5.4 


53.1 


101.7 


160.1 




2 





53.7 


65.9 


23.6 




11.4 






15.4 






2 


1 


34.5 


38.4 


4.6 


5.7 


6.8 


5.7 


9.9 


7.0 


13.4 




2 


2 


21.5 


22.2 


12.3 


2.7 


2.3 


3.6 


2.2 


15.2 


14.2 




3 









17.1 




9.8 




78.5 


182.7 






3 


1 


2.3 


2.0 


10.4 


7.2 


3.5 


5.3 


49.8 


81.1 


33.3 




3 


2 


1.2 


0.6 


3.0 


2.1 


1.3 


1.7 


13.3 


26.7 


16.3 




3 


3 


0.5 


0.9 


1.2 


1.3 


0.8 


0.9 


3.3 


8.9 


10.2 



Table 3: Total variation (maximum minus minimum) of the reduced response functions as 
a percentage of the variation of the unpolarized longitudinal. The responses are computed in 
DWIA for q = 500 Me V/c. 
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w[MeV] 


J 


M 


W L 


W T 


W TL 


W TL 


W TT 


W TT 


W T ' 


W TL ' 


w TL ' 


180 








100.0 


162.6 


30.4 




3.6 






22.9 






1 









1.4 




7.3 




94.0 


181.7 






1 


1 


5.9 


31.0 


12.0 


20.1 


12.0 


6.2 


105.2 


177.9 


230.6 




2 





80.9 


112.5 


10.7 




10.1 






38.0 






2 


1 


49.8 


80.0 


9.5 


13.0 


7.3 


6.4 


8.5 


23.6 


12.8 




2 


2 


14.5 


21.2 


8.0 


3.8 


3.7 


3.6 


3.3 


2.3 


1.2 




3 









37.0 




19.3 




246.0 


271.1 






3 


1 


1.2 


6.7 


23.3 


18.3 


8.5 


11.0 


76.4 


187.3 


76.1 




3 


2 


0.4 


1.5 


3.7 


3.4 


3.4 


5.1 


18.9 


29.4 


27.1 




3 


3 


0.2 


0.1 


0.6 


0.7 


0.7 


0.8 


4.9 


4.1 


7.5 


241 








100.0 


177.0 


45.7 




6.2 






80.4 






1 









14.2 




12.8 




164.8 


42.0 






1 


1 


43.7 


80.1 


10.4 


32.5 


5.8 


3.3 


29.7 


117.4 


212.9 




2 





53.2 


99.5 


22.8 




3.9 






5.7 






2 


1 


5.0 


15.6 


7.8 


0.5 


4.4 


6.2 


20.6 


11.9 


6.5 




2 


2 


23.5 


41.0 


12.0 


1.7 


2.0 


2.0 


5.5 


28.9 


19.3 




3 









3.4 




10.0 




89.9 


43.2 






3 


1 


2.1 


5.2 


11.4 


2.5 


0.5 


0.4 


26.0 


72.3 


28.4 




3 


2 


0.9 


2.0 


0.5 


1.5 


1.8 


1.9 


15.8 


7.8 


4.5 




3 


3 


0.7 


2.2 


0.7 


0.9 


0.5 


0.5 


2.8 


14.9 


15.0 


300 








100.0 


188.4 


40.1 




5.2 






91.0 






1 









13.5 




1.9 




174.5 


170.8 






1 


1 


52.4 


86.6 


10.4 


38.5 


6.3 


3.0 


127.4 


206.8 


223.9 




2 





80.7 


177.1 


9.4 




15.9 






59.4 






2 


1 


47.3 


92.4 


7.7 


10.8 


10.1 


8.3 


19.8 


33.4 


40.8 




2 


2 


16.4 


31.2 


11.1 


4.1 


5.3 


5.9 


11.2 


26.3 


16.2 




3 









16.3 




4.8 




294.0 


335.3 






3 


1 


6.0 


9.7 


8.9 


4.6 


1.3 


3.0 


118.5 


195.1 


93.4 




3 


2 


2.4 


5.1 


1.8 


0.9 


0.8 


0.6 


25.1 


38.2 


30.0 




3 


3 


0.6 


1.7 


0.7 


1.1 


0.5 


0.5 


6.1 


8.1 


10.8 



Table 4: Total variation (maximum minus minimum) of the reduced response functions as 
a percentage of the variation of the unpolarized longitudinal. The responses are computed in 
D WIA for q = 700 Me V/c. 
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Figure 1 : Response function Wqq + ^ shown as a function of the missing momentum p for the 
nine kinematical conditions specified in Table 1. Solid lines: computed with the complete 
optical potential; dashed lines: computed with the spin-orbit potential equal to zero; dot- 
dashed lines: computed with the total potential and with Perey parameter (3 = 0.85 in the 
wave function; dotted curves: PWIA. 



Figure 2: The same as Fig. 1, but for the vector longitudinal response W^} ^ 



Figure 3: The same as Fig. 1, but for the quadrupolc longitudinal response W^q^ 



Figure 4: The same as Fig. 1, but for the unpolarized transverse response W 00 



Figure 5: The same as Fig. 1, but for the unpolarized transverse- longitudinal response 
<o L(+) 



Figure 6: The same as Fig. 1, but for the vector TL response W 10 v 1 



T L(-\-) 

Figure 7: The same as Fig. 1, but for the quadrupole TL response W 20 
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Figure 8: The same as Fig. 1, but for the vector T' response W 10 y 1 



Figure 9: The same as Fig. 1, but for the vector T' response W n 



t' ( — ) 

Figure 10: The same as Fig. 1, but for the octupole T' response W 30 ■ 1 



Figure 11: The same as Fig. 1, but for the octupole T' response W 31 



Figure 12: The same as Fig. 1, but for the fifth response function W 00 K 1 



Figure 13: The same as Fig. 1, but for the vector TV response W w 



Figure 14: The same as Fig. 1, but for the quadrupole TV response ^ 



Figure 15: The same as Fig. 1, but for the octupole TV response W 30 



Figure 16: Response function W 0Q computed in DWIA with the complete optical poten- 
tial, shown as a function of the missing momentum p, for the nine kinematical conditions 
specified in Table 1. Solid lines: total response including all the pieces of the electromag- 
netic current; dashed lines: only including charge plus magnetization (zeroth-order) terms; 
dot-dashed lines: including, in addition, the convection current. 



Figure 17: The same as Fig. 16, but for the vector TL response W^f '. 
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